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1  Executive  Summary 

Metamaterials  such  as  left-handed  materials  (LHMs),  where  both  permittivity  and  permeability 
are  negative,  have  attracted  great  attention  due  to  their  novel  and  unique  electromagnetic  proper¬ 
ties.  The  metamaterials  have  potential  applications  such  as  radio  frequency  and  microwave  device 
miniaturization.  Because  excited  complex  surface  modes  have  slowly  exponential  decaying  fields  on 
the  air-slab  interface,  the  directivity  for  dipoles  over  a  grounded  metamaterial  slab  is  much  higher 
than  directivity  with  conventional  double  positive  (DPS)  materials.  This  structure  is  very  useful 
to  study  slot  antennas  and  leaky-wave  antennas. 

We  propose  to  investigate  leaky  waves  and  their  excitation  with  multilayered  metamaterials. 
The  objective  of  this  project  is  to  enhance  antenna  performance  with  excitations  of  leaky  waves. 
First,  we  investigate  a  grounded  dielectric  slab  with  double  negative  (DNG)  materials.  Dramatically 
different  dispersion  curves  of  evanescent  surface  modes  (electromagnetic  fields  that  exponentially 
decay  both  in  air  and  inside  the  slab)  are  observed,  indicating  that  they  are  highly  dependent  on 
the  medium  parameters.  As  the  counterpart  of  the  improper  complex  leaky  modes  in  a  slab  with 
a  double  positive  (DPS)  medium,  the  complex  modes  in  a  slab  with  a  DNG  medium  are  proven 
to  be  exclusively  proper.  They  have  exponentially  decaying  fields  in  the  air  region  and  are  termed 
complex  surface  modes.  There  are  an  infinite  number  of  complex  surface  modes  and  they  cannot 
be  suppressed.  Further  investigation  of  the  Poynting  vector  shows  that  those  complex  modes  do 
not  carry  away  power  in  both  transverse  and  longitudinal  directions.  However,  the  complex  surface 
modes  play  an  important  role  in  the  far  field  radiating  from  the  interface  between  the  air  region 
and  the  grounded  slab. 

Secondly,  we  read  Tamil'  and  Kou’s  paper  and  other  papers  on  leaky  waves  in  asymmetric  slabs. 
By  ’Asymmetric’  it  means  that  the  material  above  and  below  the  slab  can  have  different  dielectric 
constant,  which  makes  the  discussion  more  general.  Tarnir  and  Kou  claim  that  there  are  eight 
different  leaky-wave  fields  guided  by  an  asymmetric  layer  configuration,  but  only  four  of  the  eight 
can  be  excited  in  conventional  DPS  materials.  Those  eight  leaky  waves  can  be  divided  into  four 
forward  and  four  backward  leaky  waves,  or  four  co-  and  four  contra-leaky  waves.  It  is  found  that 
for  a  lossless  DPS  grounded  slab,  there  are  only  two  improper  complex  modes  (leaky  waves).  If  the 
slab  is  lossy,  the  conventional  surface  wave  modes  become  complex  modes,  which  are  called  leaky 
waves  by  Tamir  and  Kou.  A  plasma  slab  will  add  another  type  of  proper  complex  modes.  For  a 
grounded  slab  with  DNG  materials,  only  two  proper  regular  complex  modes  exist  for  a  lossless  slab. 
If  the  slab  is  lossy,  the  proper  surface  wave  modes  may  become  proper  complex  modes  or  improper 
complex  modes. 

Thirdly,  we  investigated  the  Sommerfeld  integral  path  (SIP)  for  dielectric  slab  with  double 


1 


negative  metamaterials.  It  is  found  that  layered  double  negative  materials  have  different  SIPs  as 
compared  to  conventional  layered  double  positive  media.  In  a  layered  medium,  the  Green’s  function 
is  written  as  a  Sommerfeld-type  integral,  whose  integrand  has  poles  on  the  hp-plane  which  corre¬ 
spond  to  the  eigenmode  solutions  of  the  same  structure  with  the  source  removed.  After  the  proper 
integral  path  has  been  established,  there  are  two  ways  to  represent  the  total  hied.  The  first  one  is 
often  to  deform  the  path  to  a  branch  cut  (BC)  with  poles  captured  by  the  branch  cut,  which  leads 
to  a  spectral  representation.  The  second  way  is  to  deform  the  path  to  a  steepest  descent  path  (SDP) 
with  poles  captured,  which  lead  to  a  non-spectral  representation  for  calculating  far  held  (from  the 
integral  along  the  SDP  directly).  When  the  medium  is  lossless,  the  Sommerfeld  integral  blows  up 
at  real  surface  poles  along  the  real  axis.  A  common  way  to  address  this  difficulty  is  to  assume  that 
the  medium  is  lossy  and  then  monitor  how  pole  loci  change  as  the  loss  becomes  smaller  and  smaller. 
From  our  previous  results,  the  Sommerfeld  integral  path  is  chosen  as  a  half  circle  above  or  below  the 
poles  depending  on  the  kind  of  poles,  the  hrst  kind  or  second  kind.  Physically  the  hrst  kind  of  poles 
corresponds  to  forward  surface  modes  (the  power  how  is  in  the  same  direction  as  the  wave  propa¬ 
gation),  while  the  second  kind  of  poles  corresponds  to  backward  surface  modes  (the  power  how  is 
in  the  opposite  direction  as  the  wave  propagation).  We’ll  talk  about  the  poles  in  details  in  Section  3. 

Finally,  we  investigated  the  radiation  intensity  and  directivity  of  electric/magnetic  dipoles  over 
a  grounded  slab  to  find  optimal  values  of  material  parameters  and  thickness  for  a  given  frequency. 
We  considered  a  horizontal  magnetic  dipole  (HMD)  placed  on  a  PEC  and  covered  by  a  layer  of 
DNG  material  with  thickness  d  and  constitutive  parameters  er  and  /ir .  The  material  is  homoge¬ 
nous,  isotropic,  non-dispersive.  We  optimize  the  directivity  over  the  thickness  and  the  constitutive 
parameters.  It  is  found  that  for  a  slab  with  a  hxed  electric  size  (proportional  to  the  product  of 
the  thickness  and  frequency),  the  directivity  reaches  the  maximum  when  the  product  of  er  and 
fir  is  a  certain  value,  which  is  determined  by  the  electric  size  only.  The  maximum  directivity  is 
proportional  to  er  and  {d/ Ao)8^3.  We  also  found  the  fitting  formula  for  the  maximum  directivity 
and  /ir  as  functions  of  the  electric  size  and  er. 

This  final  project  report  is  organized  as  follows:  The  second  section  gives  a  brief  introduction  to 
this  field.  Section  3  briefly  reviews  the  mode  properties,  dispersion  diagrams,  and  power  flows  and 
plots  pole  loci  versus  frequency  and  media  loss.  Section  4  discusses  the  proposed  new  Sommerfeld 
integral  path  and  addresses  two  numerical  issues  associated  with  Sommerfeld  integral  that  are 
generally  overlooked  by  the  literature.  Section  5  gives  several  numerical  examples  that  demonstrate 
the  necessity  of  the  new  SIP  for  DNG  media,  and  the  superior  radiation  patterns  of  a  simple  dipole 
antenna  filled  with  DNG  media.  We  also  optimize  structure  for  maximum  directivity  of  HMD  at 
the  broadside.  Finally,  Section  6  concludes  the  report  with  a  summary,  publication  list,  significant 
events,  and  proposed  future  work. 
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2  Introduction 


Due  to  the  progress  in  fabrication  technology  and  nanotechnology,  there  is  a  renewed  interest  in 
developing  composite  materials  that  mimic  known  material  responses  or  that  qualitatively  have 
new,  physically  realizable  response  functions  that  do  not  occur  in  nature  [1].  Examples  of  engi¬ 
neered  materials  include  left-hand  materials  (LHMs),  negative  refraction  materials,  electromagnetic 
band-gap  (EBG)  structured  materials,  and  complex  impedance  surfaces.  Those  new  metamaterials 
often  are  generated  by  artificially  fabricated  inhomogeneities  embedded  in  host  media. 

The  metamaterials  have  potential  applications  such  as  radio  frequency  and  microwave  device 
miniaturization.  Microstrip  antennas  are  easy  to  make  and  have  the  advantages  of  small  size,  light 
weight,  and  low  cost.  Surface  waves  always  exist  on  grounded  conventional  dielectric  slabs,  carry 
away  energy,  and  reduce  antenna  directivity.  Leaky  waves  are  a  type  of  complex  surface  waves 
and  have  slowly  exponential  decaying  fields  along  the  air-slab  interface.  The  leaky  waves  make  the 
fields  on  the  interface  equivalent  to  an  antenna  with  large  aperture,  which  leads  to  the  directivity 
for  dipoles  over  a  grounded  metamaterial  slab  is  much  higher  than  directivity  with  conventional 
double  positive  (DPS)  materials.  This  structure  is  very  useful  to  study  slot  antennas  and  leaky- 
wave  antennas. 

Pendry  et  al.  proposed  forming  LffMs  by  introducing  a  periodic  array  of  split  ring  resonators 
and  continuous  wires  [2],  LHMs  have  been  demonstrated  experimentally  by  Smith  et  al.  and 
Shelby  et  al.  [3,4],  Double  negative  (DNG)  metamaterials  have  attracted  intensive  interest  in  the 
last  few  years  for  their  exotic  properties,  such  as  negative  refraction,  reversed  Doppler  effect,  re¬ 
versed  Vavilov-Cerenkov  effect,  and  the  possibility  of  making  perfect  lens  [5,6,2,7,3,8,9,10,11, 
12,13,14,15].  Many  researchers  have  shown  great  interest  in  antenna  applications  of  metamate¬ 
rials  [1,16,17,18,19,20].  Ziolkowski  and  Kipple  [16]  showed  analytically  that  the  LHMs  increase 
radiation  efficiency  of  small  simple  antennas.  An  electric  or  magnetic  dipole  radiating  in  a  grounded 
dielectric  slab  with  ideal  LHMs  has  been  studied  numerically  and  the  proper  complex  modes  (leaky 
waves)  result  in  a  radiation  pattern  like  a  pencil  beam  at  broadside  [21,22], 

Among  the  many  applications  and  novel  structures  proposed  for  double  negative  metamateri¬ 
als,  the  grounded  dielectric  slab  with  a  DNG  medium  is  the  most  fundamental  structure  because 
it  not  only  is  the  theoretical  foundations  for  other  more  complicated  structures  but  also  has  many 
potential  applications  such  as  transmission  line  in  real  systems.  [23]  and  [24]  found  that  there  are 
special  regions  for  transverse  magnetic  (TM)  modes  where  two  different  propagation  constants  ex¬ 
ist.  [25,17,26,27]  studied  evanescent  surface  modes,  which  do  not  exist  in  a  double  positive  (DPS) 
medium.  The  grounded  dielectric  slab  also  supports  complex  surface  modes  even  for  lossless  me¬ 
dia  [28,18,27].  Baccarclli  and  his  colleagues  proposed  the  concept  of  surface  wave  suppression, 
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which  ensures  the  absence  of  both  ordinary  and  evanescent  surface  modes.  This  is  very  attractive 
in  view  of  taking  the  DNG  medium  as  a  potential  substrate  candidate  to  reduce  edge  diffraction 
effects  and  enhance  radiation  efficiency  for  microstrip  antennas  [29,21].  Later,  the  necessary  and 
sufficient  conditions  for  surface  wave  suppressions  were  derived  in  [30,27]. 

In  [18],  radiations  from  a  magnetic  dipole  on  grounded  dielectric  slab  with  DNG  medium  are 
studied.  In  analyzing  such  problems,  the  spectral  domain  approach  is  used  and  consists  of  two 
steps:  first,  the  electromagnetic  fields  generated  by  a  dipole  are  solved  in  the  spectral  domain  via 
a  Fourier  transform;  the  spectral  domain  solution  is  then  transformed  back  to  the  space  domain. 
The  second  step  involves  a  Sommerfcld-type  integral  (integral  from  zero  to  infinity).  For  a  lossless 
medium  there  are  a  branch-cut  issue  and  surface  mode  singularity  blow-up  problem  [31].  A  com¬ 
mon  way  to  remedy  the  aforementioned  difficulties  is  to  firstly  assume  that  dielectric  medium  is 
lossy  and  then  let  the  loss  decrease  to  zero  [32].  With  this  method,  the  Sommerfeld  integral  path 
(SIP)  for  conventional  lossless  medium  is  well-defined  [31]. 

However,  the  grounded  dielectric  slab  with  a  DNG  medium  has  a  mode  spectrum  that  is  much 
more  complicated  than  that  of  a  conventional  double  positive  medium  [27] .  Its  complete  eigenmode 
spectrum  consists  of  evanescent  surface  modes,  ordinary  surface  modes,  complex  surface  modes,  and 
continuous  radiation  modes.  By  contrast,  the  complete  mode  spectrum  of  a  DPS  medium  consists 
of  only  ordinary  surface  modes  and  radiation  modes.  In  addition  to  the  mode  expansion,  the 
ordinary  surface  modes  of  a  DNG  medium  bifurcate  near  their  cutoff  frequencies.  Both  of  these  two 
phenomena  contribute  to  the  novel  properties  of  the  DNG  medium.  On  the  other  hand, they  make 
the  problems  related  to  DNG  medium  more  complicated.  One  of  the  most  significant  ramifications 
is  that  a  DNG  medium  requires  a  different  Sommerfeld  integral  path,  which  will  be  elaborated  on 
in  this  report  in  great  detail. 
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Mode  Properties  and  Surface  Pole  Loci 


3.1  Eigen  Value  Problem 

The  structure  we  investigate  is  a  grounded  dielectric  slab  of  thickness  d  (see  Figure  1).  Region  I 
is  air  with  er\  and  fiT \  and  Region  II  is  a  DNG  medium  er2  and  fir2 .  Inherently,  the  passive  DNG 
media  must  be  dispersive  [9,33].  However,  as  is  common  in  the  literature,  most  of  the  DNG  medium 
considered  in  this  paper  is  assumed  to  be  non-dispersive  for  simplicity  purpose.  This  assumption 
is  found  acceptable  because  a  small  perturbation  of  er  and  /ir  satisfies  its  dispersiveness. 


Crlj  Rrl  ®  ^  V 

X 

z  =  d 


z  =  0 


Figure  1:  A  grounded  dielectric  slab  with  lossless,  isotropic  and  homogeneous  DNG  medium. 


3.1.1  Decoupling  of  TE  and  TM  Modes 

For  homogeneous  isotropic  medium,  the  two  curl  equations  in  source  free  region 


(1) 

(2) 


V  x  E  =  —  juif-tJi 

V  x  H  =  jive E 
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are  expanded  into  six  scalar  equations 


dEy 

dz 


+  lyEz  — 


^)fyEx 


dHy 

dz 


d  Ex  _ 
dz 

T 


^yHx 
dHx  _ 
dz 

In  deriving  the  above  equations,  d/dx  —  0  and  d/dy 


ju>nHx 

(3a) 

ju/iHz 

(3b) 

jujyHy 

(3c) 

-jcveEx 

(3d) 

jueEz 

(3e) 

jcoeEy 

(3f) 

=  — t /y  are  used. 


Equations  (3a),  (3e),  and  (3f)  are  functions  of  Ey,  Ez,  and  Hx  only.  And  (3b),  (3c),  and 
(3d)  are  functions  of  Hy,  Hz,  and  Ex  only.  This  observation  implies  that  under  the  condition  of 
d/dx  =  0,  the  electromagnetic  waves  are  decoupled  into  independent  TE  modes  ( Hy ,  Hz,  Ex )  and 
TM  modes  ( Ey ,  EZJ  Hx).  Notice  that  in  this  chapter,  TE  and  TM  modes  are  referred  to  horizontal 
propagation  direction  (^/-direction  in  Figure  1).  This  classification  is  also  applicable  to  vertical 
direction  (^-direction  in  Figure  1)  as  will  be  seen  in  later  chapters. 


3.1.2  Homogeneous  Helmholz  Equation 

Substituting  (3b)  and  (3c)  into  (3d),  we  arrive  at  1-D  homogeneous  Helmholz  equation  for  TE 
modes  propagating  along  ^/-direction 


where  kf  =  u2eifii,  i  =  1,2,  and 


d2Ex 
dz 2 


+  (ki  + 


0 


(4) 


7  y  =  &  +  jfd 

With  the  same  procedure,  we  have  1-D  homogeneous  Helmholz  equation  for  TM  modes 


(5) 


d2Hx 

dz2 


+  (k2  +  ^ 


0 


(6) 
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3.1.3  Boundary  Conditions 

Equations  (4)  and  (6)  are  1-D  differential  equations  valid  for  each  source  free  region.  Boundary 
conditions  must  be  applied  to  solve  the  problem. 


At  the  interface  of  Region  1  and  Region  2,  tangential  electric  and  magnetic  fields  are  continuous. 
So  for  TE  modes 


|  z=d 
Hyl  |  z=d 

For  TM  modes 

Aj/l  \z=d 
Hx  1  |  z=d 


Ex2  |  z=d  (7a) 

Hy2  I  z=d  (7b) 

77^2  |  z=d  (8a) 

Hx 2  | z=d  (8b) 


In  Region  1,  proper  boundary  condition  requires  vanishing  field  at  infinity 


A 

I  z=oo 


0 


(9) 


or  radiation  condition  [34] 

lim  r 

r— >oo 

where  A  =  ExA  for  TE  modes  and  A  =  HxA  for  TM  modes.  Notice  that  (9)  is  a  special  case  of  (10). 
Solutions  do  not  satisfy  (10)  are  non-physical.  However,  they  are  still  of  importance  as  we  will  see 
later. 


<9A  ' 

vp  +  jkA 
or 


=  0 


(10) 


In  Region  2,  there  is  boundary  conditions  of  PEC  at  z  =  0.  For  TE  modes 


For  TM  modes 


2  =  0 


0 


dHx  2 


dz 


2=0 


0 


(11) 

(12) 


3.1.4  Discrete  TE  Modes 

The  general  solutions  of  (4)  in  region  1  and  region  2  are 

Exl  =  Ae~^z-d)  +  Be^l{z~d)  (13) 

Ex2  =  Ce~^2Z  +  De^2Z  (14) 
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where  y2i  =  —klerin,ri  —  7^  and  i  =  1,2  is  the  layer  index. 

To  satisfy  boundary  condition  of  (11),  we  have  C  =  —D  or 

Ex2  =  C'sinh(7*2  z)  (15) 

Surface  waves  have  zero  held  at  infinity.  Thus  in  (13),  one  term  must  be  removed.  The 
choice  is  arbitrary  since  we  have  two  different  choices  for  7^1  too.  For  example,  one  can  choose 
Ex  1  =  Beyzl('z-d'>  and  7^  =  —  — fcgeri//ri  —  7^.  This  is  the  same  as  the  choice  of  Exl  =  Ae~lzl('z~d') 

and  7^1  =  yj— kleri/iri  —  7 y.  To  be  consistent  with  literature,  we  use  Exl  =  yfe-721^-^  and  assume 
7^1  has  positive  real  part. 

With  the  above  observation,  Ex  and  Hy  in  each  region  are 


Ex  1  =  Ae~^l{z~d)  (16a) 

Ex2  =  Csmh(jz2z)  (16b) 

Hy\  = - —  =  J^_Ae-^l{z~d)  (16c) 

JUVi  oz  juni 

Hy2  =  — — —  =  — ^-C  cosh  (7z2  z)  (16d) 

ju/d2  dz  ju/d2 

Using  boundary  condition  of  (7),  we  have 


1  -  sinh(722d)  1  [  A  1  _ 

251  ^  cosh  (7 z2d)  C  ' 

M  Mr  2  V  IZZ  >  J  L  J 

Non-trivial  solutions  require  the  determinant  of  the  matrix  to  be  zero.  Therefore  we  arrive  at  eigen 
equation  for  TE  modes 

—  sinh(7z2d)  +  —  cosh(7z2d)  =  0  (17) 

/iri  /ir2 

where  qz2  7^  0.  To  implicitly  eliminate  the  fake  solution  of  yz2  =  0,  (17)  is  written  as 

—  sinh(722d)  +  —  cosh(7z2d)  =  0  (18) 

7;j2  Iz  1 


3.1.5  Discrete  TM  Modes 

With  the  same  procedure,  eigen  equation  for  TM  modes  is 

—  sinh(7*2d)  +  —  cosh(722d)  —  0 


(19) 


(20a) 

(20b) 


(20d) 


3.1.6  Field  Components  of  Discrete  TE  and  TM  Modes 

The  field  components  of  discrete  TE  modes  are  (ryZ2  ^  0) 

Exi(y,z)  =  —  sinh(722 d)e~'yzl('z~d'>e~lyV  (20a) 

Ex2{y,z)  =  -sinh(7z2^)e-7^  (20b) 

Hyl(y,z)  =  -.1  j_  sinh(7 z2d)e-^z-^e-™  (20c) 

juyi  dz  ju/i i 

Hy2(y,  z )  =  — r~  =  — ~  cosh(722^)e-7^  (20d) 

juy2  dz  juy2 

Hzl{y ,  z)  =  ~^Exl  =  sinh(7 z2d)e~^-<D  e-ivy  (20e) 

juyi  juyi 

Hz2(y ,  z)  =  — t —Ex2  =  -EL-  sinh(722^)e_7^  (20f) 

juy2  juy  2 

The  field  components  of  discrete  TM  modes  are 

Hxl(y,z)  =  -cosh(722d)e_7zl(z_,i)e_7^  (21a) 

Hx2(y,  z)  =  -  cosh(jz2z)e~lyy  (21b) 

Eyi(y,  z)  =  — —  =  -EL  cosi^7z2^e-7zi(z-d)e-7yy  (2ic) 

jcuei  dz  ]uex 

Ey2(y,  z )  =  — —  =  -^LL  sinh(7z2^)e_7y?/  (21d) 

]ue2  dz  jue2 

Ezl(y ,  z)  =  ~^—Exl  =  -^cosh(lz2d)e-^-d)e-7yy  (21e) 

juei  jue  i 

Ez2{y,  z)  =  ——  Ex2  =  — cosh(722^)e-7^  (21f) 

jue2  jue2 

3.1.7  Graphical  Method 

For  ordinary  real  surface  modes,  7,1  =  azl  and  722  =  j/3z 2,  where  both  a2l  and  fdz2  are  real  numbers. 
Therefore  (18)  and  (19)  are  rewritten  as  follows 

^— ■( /3z2d, )  cot (/3z2d)  =  —az\d  for  TE  mode  (22) 

t^r2 

—  (/3z2d)t&ii(l3z2d)  =  az\d  for  TM  mode  (23) 

Cr-2 

Furthermore,  using  relations  7^  =  —k^erifiri  —  L  and  i  —  1,2  and  eliminating  7^  lead  to 


(a2id)2  +  (fiz2d)2  —  (kod)2  (er2yr2  —  er\yri) 
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The  evanescent  surface  modes  have  exponentially  decaying  electromagnetic  fields  in  both  air 
and  dielectric  layer.  For  evanescent  surface  modes,  7-1  =  az\  and  7^2  =  ay 2,  where  both  az\  and 
az 2  are  real  numbers.  The  eigen  equations  (18)  and  (19)  are  rewritten  as 


(a22<7)  coth(ay2d )  =  —az\d 

P-r2 

for  TE  mode 

(25) 

—  (az2d)  tanh(ay2G0  =  —azld 
er2 

for  TM  mode 

(26) 

The  constitutive  equations  are  now 

(a2ld)  (a2 2^)"  (A^od)  (er-2/^r2  ^r\  /A  1 )  (27) 

The  graphical  representations  [35,36]  of  the  above  equations  are  shown  in  Figure  2.  In  fact  the 
figure  combines  two  figures.  The  right-half  is  for  the  ordinary  modes  governed  by  Equations  (22)  to 
(24),  and  the  left  half  is  for  the  evanescent  modes  governed  by  Equations  (25)  to  (27).  The  index 
notation  used  here  follows  [35]:  only  TE  even  modes  (but  with  odd  indices,  TEx,  TE3,  . . . )  and  TM 
odd  modes  (but  with  even  indices,  TM0,  TM2,  . . . )  can  exist  in  a  PEC  grounded  dielectric  slab. 
Notice  that  in  the  first  and  second  quadrants,  az\  is  positive  and  the  fields  exponentially  decay  in 
the  air  region  (proper);  in  the  third  and  fourth  quadrants,  az\  is  negative  and  the  fields  exponen¬ 
tially  increase  in  the  air  region  (improper).  The  x-axis  is  divided  into  two  segments.  The  right 
half  is  for  /3z2 d,  whose  fields  inside  the  dielectric  layer  are  sine/cosine  standing  waves  (ordinary); 
while  the  left  half  is  for  ay2d,  whose  fields  inside  the  dielectric  layer  are  exponentially  distributed 
(evanescent).  The  intersection  in  the  second  quadrant  represents  the  proper  evanescent  surface 
mode  which  does  not  exist  for  a  DPS  medium.  This  can  be  seen  from  (22),  (23),  (25),  and  (26): 
the  difference  for  DPS  and  DNG  media  is  a  minus  sign  on  the  left  sides  of  these  equations,  which  im¬ 
plies  that  one  can  get  results  for  a  DPS  medium  by  mirroring  Figure  2  along  its  horizontal  axis.  By 
doing  this,  the  intersection  in  the  second  quadrant  will  be  in  the  third  quadrant,  and  thus  improper. 

Another  important  fact  that  must  be  recognized  from  Figure  2  is  that  the  ordinary  surface 
wave  solutions  for  a  DNG  medium  are  no  longer  monotonic.  It  is  clear  from  the  sub-figure  in 
the  left  corner  of  Figure  2(a)  that  there  are  two  intersections  as  the  radius  of  the  dashed  circle 
decreases,  which  corresponds  to  a  decrease  in  electric  size,  i.e.,  k0d.  Once  the  circle  has  only  one 
tangential  point  with  the  solid  line,  further  decreasing  k^d  will  cause  this  mode  to  be  cutoff.  Similar 
occurrence  happens  to  TM  modes  in  Figure  2(b)  in  a  more  obvious  way.  These  two  possible  modes 
have  two  different  power  distributions.  One  has  more  power  flowing  in  the  air  region  than  that  in 
the  dielectric  region,  making  the  total  power  flow  in  the  same  direction  as  the  phase  velocity.  The 
other  is  in  the  opposite  way  and  displays  the  backward  property.  More  details  on  the  Poynting 
vectors  are  addressed  in  later  sections. 
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47r  37r  27r  7r  0  7r  2ir  37 r  47t 


(a)  TE  modes 


Figure  2:  Graphical  solutions  for  TE  and  TM  modes.  Solid  lines  in  the  first  and  fourth  quadrants 
represent  (22)  or  (23);  solid  lines  in  the  second  quadrant  represent  (25)  or  (26);  dashed  line  in  the 
first  and  fourth  quadrants  represents  (24);  dashed  line  in  the  second  and  third  quadrants  represents 
(27).  The  medium  parameters  are:  er\  =  1,  /ir i  =  1,  er2  =  —2.5,  /rr2  =  —0.5. 


11 


3.2  Mode  Properties 

This  section  states  different  mode  properties,  such  as  evanescent  surface  modes,  complex  surface 
modes,  Poynting  vectors,  etc.  We  define  some  terms  here  first. 

Discrete  or  continuous  modes:  Discrete  modes  refer  to  discrete  solutions  of  eigen  equations 
derived  in  last  section  and  continuous  modes  are  the  waves  satisfying  the  radiation  boundary  condi¬ 
tion  and  having  continuous  spectrum,  which  are  the  contribution  for  the  integral  along  the  branch 
cuts. 

Proper  or  improper  modes:  Proper  modes  exist  physically  and  the  fields  of  the  improper 
modes  increase  exponentially  to  infinity  at  the  infinity.  So  the  improper  modes  cannot  exist  phys¬ 
ically  and  are  not  part  of  complete  spectrum.  For  the  improper  modes,  az\  =  Re^i]  <  0.  For 
the  proper  modes,  Re  [771]  >  0.  The  modes  with  Re[7zl]  >  0  are  also  called  surface  wave  modes 
because  the  fields  decay  exponentially  away  from  the  interface  in  the  top  region. 

Real  or  complex  modes:  For  a  lossless  slab,  the  waves  with  a  real  propagation  constant  (7^ 
or  7 p)  are  called  as  real  modes,  otherwise,  they  are  called  as  complex  modes.  In  general,  all  modes 
are  complex  for  a  lossy  slab. 

Ordinary  or  evanescent  modes:  Two  types  of  real  modes  may  be  supported  by  a  lossless 
slab:  ordinary  modes  (standing  waves  in  the  dielectric  region,  072  =  Re[7Z2]  =  0)  and  evanescent 
modes  (the  fields  change  exponentially  in  the  dielectric  region,  f5z 2  =  Im[7z2]  =  0). 


3.2.1  Evanescent  Surface  Wave  Modes 

The  normalized  effective  dielectric  constant  ee//  =  ((3/ko)2  for  evanescent  surface  wave  modes  is 
larger  than  both  er\nr\  and  er2/ir2-  Therefore  yZ2  =  \J — —  7y  =  &oa Ae//  —  er2/ir2  is  a  pure 
real  number.  The  electromagnetic  fields  inside  the  dielectric  layer  are  no  longer  sine/cosine  standing 
waves,  but  of  the  form  of  A  sinh(aZ2^)-  Figure  3  shows  a  typical  held  configuration  for  evanescent 
surface  wave  modes.  Notice  that  the  held  extends  to  the  air  region  far  away  and  decay  very  slowly. 

It  is  found  that  the  dispersion  curves  for  evanescent  surface  wave  modes  are  complicated  and 
highly  dependent  on  medium  parameters.  Figure  4  shows  two  dispersion  diagrams  for  TEi  mode 
with  different  medium  parameters.  The  solid  line  in  Figure  4  is  for  proper  mode  while  the  dotted 
line  is  for  improper  mode.  The  dashed  lines  in  both  figures  depict  the  value  of  y/er2fir2-  They  are 
the  watersheds  by  which  one  can  tell  evanescent  surface  wave  modes  from  ordinary  ones. 
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Figure  3:  Demonstration  of  field  distribution  for  a  TE  evanescent  surface  wave  mode.  Medium 
and  structure  parameters  are:  er\  =  1,  fir i  =  1,  er2  =  —2.5,  /ir2  =  —2.5,  k^d  =  0.1  rad.  Solved 
propagation  constant:  =  j4. 27819. 


In  Figure  4(a),  evanescent  surface  wave  mode  has  low  cutoff  point.  As  the  electric  size  in¬ 
creases,  the  ordinary  surface  wave  mode  becomes  an  evanescent  surface  wave  mode  and  its  effective 
dielectric  constant,  ee//,  keeps  increasing.  In  Figure  4(b),  however,  the  situation  is  reversed.  The 
evanescent  surface  mode  has  a  high  cutoff  point  above  which  it  becomes  the  ordinary  surface  wave 
mode.  When  k0d  is  small,  the  evanescent  surface  wave  mode  has  an  extremely  large  eeff  that 
decreases  rapidly  as  the  electric  size  increases.  One  can  refer  to  the  sub-figures  of  Figure  4  to  check 
its  validations,  which  are  graphic  solutions  shown  in  earlier  section. 


3.2.2  Complex  Surface  Modes  and  Its  Poynting  Vectors 

It  is  well  known  that  the  complete  proper  mode  spectrum  of  a  lossless  DPS  dielectric  slab  include 
discrete  surface  wave  modes  and  continuous  radiation  modes,  both  of  which  are  real  surface  wave 
modes  [34,36].  With  a  DNG  medium,  it  is  proved  that  all  complex  roots  of  the  eigen  equations 
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Figure  4:  Two  possible  dispersion  curves  for  TE  proper  surface  modes  (solid  lines)  and  TE  improper 
leaky  modes  (dotted  lines).  The  dashed  line,  representing  ^/er2/rr2,  is  the  watershed  for  evanescent 
surface  mode  and  ordinary  surface  modes. 
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Figure  5:  Demonstration  of  field  distribution  for  a  TE  complex  surface  wave  mode.  Medium  and 
structure  parameters  are:  eri  =  1,  Hr i  —  1,  er2  =  —2.5,  /xr 2  =  —2.5,  k0d  =  1.5  rad.  Solved 
propagation  constant:  7y//co  =  0.8261  +  j  1. 1253. 


are  on  the  top  Riemann  sheet.  These  solutions,  termed  complex  surface  waves,  form  another  set 
of  proper  modes.  Unlike  real  surface  wave  modes,  complex  surface  wave  modes  have  high  cutoff 
electric  size,  which  means  that  they  exist  only  when  k^d  is  small.  Therefore,  there  are  an  infinite 
number  of  complex  surface  wave  modes  from  very  small  electric  size  all  the  way  to  large  electric 
size  as  can  be  seen  later  in  Figure  7.  Figure  5  shows  field  distribution  of  a  typical  complex  surface 
wave  mode.  The  electromagnetic  fields  decay  pretty  fast  in  the  air  region. 

The  existence  of  an  infinite  number  of  complex  surface  wave  modes  at  any  frequency  seems 
somehow  annoying  in  terms  of  suppressing  surface  wave  modes  since  they  will  carry  away  power 
and  lower  the  radiation  efficiency.  Fortunately,  a  careful  examination  on  the  Poynting  vector  elim¬ 
inates  this  possibility. 

The  Poynting  vector  is  defined  as 

S  =  x  H*  (28) 
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And  we  have  the  following  trigonometry  identities 

sinh(a  +  jb )  =  sinh  a  cos  b  +  j  cosh  a  sin  b 
cosh(a  +  jb)  =  cosh  a  cos  b  +  j  sinh  a  sin  b 

|sinh(a  +  jb)  |  =  sinh2  a  +  sin2  b  =  -  [cosh(2a)  —  cos(26)] 
|cosh(a  +  jb) |2  =  cosh2  a  —  sin2  b  =  -  [cosh(2a)  +  cos(26)] 
sinh(o  +  jb)  cosh*  (a  +  jb)  =  -  [sinh(2a)  +  j  sin(26)] 

sinh*  (a  +  jb)  cosh(a  +  jb)  =  -  [sinh(2a)  —  j  sin(26)] 

Using  the  field  components  in  Section  3.1.6,  we  have  for  TE  mode 


% 


Ql  -C/  7U  TT*  

Syl  ~  2EXlHzl- 

oTE  _  _  1  rp  TT*  _  4 y 

Ey2  -  2Ex2Hy2-  4^5 

qTE  "^771  zj*  4'z  1 

^  2  : cl^1_  4^1* 


[cosh(2a22<f)  —  cos(2 f3z2d)]  e 
[cosh(2 az2z)  -  cos(2 /3z2z)\  e~2ay 


-2az\(z-d)  e~2ay 


cosh(2az2d)  —  cos(2 fiz2d)\  e 


-2azl(z-d)  -2ay 


qTE  _  *771  TT*  _ 

bz2  ~  n^x2 ny2  - 


lz2 


2  ^'~yz  Ajuj/i* 

The  total  power  flow  in  ^-direction  for  TE  mode  is 


sinh(2 az2z)  +  j  sin(2 f3z2z)}  e  2ay 


jTE  ( 

y _ 

o-2  ay 


PTE(v)  rd 
v  =  /  S™dz  + 


dy 


1 0 


Jy  1 


% 


sinh(2  az2d)  sin(2  j3z2d) 


OLz2 


P-. 


'z2 


% 


8ju/dla 


z  1 


cosh(2  az2d)  —  cos(2  /3z2d)] 


8 jwf4 

And  for  TM  mode  we  have 

s™  =  Iev1H*x1  =  [cosh(2a22d)  +  cos(2 /3z2d)]  e-2a^z~d^ e~2ay 


4jiOCi 


Sy2M  =  -Ez2H*2  =  — [cosYi{2olz2z)  +  cos(2(3z2z)\  e  2ay 


4  jue2 
7*i 


Ql  M  _  77>  TT*  _ 

S*  -~2Ev'H*-A juei 

lz2 


cosh(2 az2d)  +  cos(2 /3z2d)\  e 


-2azi(z—d)—2ay 


qTM  _  *  rp  TT*  _ 


'  z2 


z2ELx2 


4 jue2 


sinh(2o'z2^)  +  j  sin(2 /322*)\  e  2ay 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


(39) 

(40) 

(41) 

(42) 

(43) 
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The  total  power  flow  in  ^-direction  for  TM  mode  is 


e-2ay 


S™dz  + 


S™dz 


7 y 


8juje2 


sinh(2a!22d)  sin(2  (3z2d) 


Oiz2 


Pz2 


+  0  .  - [cosh(2az2d)  +  cos(2 f3z2d)) 

8]ueiazl 


(44) 


Although  the  equations  of  the  total  power  flow  (39)  and  (44)  are  quite  complicated,  the  result 
turns  out  to  be  simply  zero  [37]  as  can  be  seen  by  the  following  physical  argument.  To  give  out 
a  clear  picture  of  the  power  flow  density,  the  Poynting  vector  is  drawn  in  Figure  6.  Dotted  lines 
are  the  boundaries  of  a  conjectured  box  with  dimension  of  4  mm  x  20  mm.  The  power  flow  has 
different  directions  in  the  air  and  inside  the  dielectric  layer.  The  average  total  power  flowing  out 
of  the  box  via  the  left  and  right  side  walls  is  equal  to  the  power  flowing  into  the  box  from  the  top 
wall.  Now  assuming  the  top  wall  is  moved  to  infinity,  there  is  no  power  flowing  into  the  box  from 
the  top  wall  since  the  complex  surface  modes  have  zero  fields  at  infinity.  According  to  the  energy 
conservation  law,  the  power  exchange  via  the  left  side  wall  Pyl  =  /0°°  Syi(z,  y)dz  must  be  equal 
to  the  power  exchange  via  the  right  side  wall  Py2  =  /0°°  Sy2(z,y)dz.  Note  that  the  propagation 
constant  along  ^-direction  is  now  a  complex  number  jy  =  a  +  j/3  where  a^O,  and  Pi  and  P2  are 
related  by  P2  =  P1e_2"^2_yib  tq  ensure  P\  =  P2,  the  only  possibility  is  P\  =  P2  =  0.  Thus  all  the 
complex  surface  modes  have  zero  power  flow  in  ^-direction.  They  do  not  transport  any  energy. 


3.3  Mode  Spectrum  of  DNG  Medium 

As  discussed  in  previous  section,  now  we  have  got  the  eigen  equations  of  transverse  electric  (TE) 
and  transverse  magnetic  (TM)  modes  for  waves  propagating  in  p-direction: 


— -  sinh(722(i)  +  cosh(722d)  =  0  TE  mode 

7*2  7*i 

(45) 

—  sinh(722d)  +  —  cosh(7z2ci)  =  0  TM  mode 

(46) 

6r2  €rl 

T/ zi  “t-  jPzi 

(47) 

where  y2l  {i  =  1,  2)  is  the  complex  propagation  constant  along  ^-direction  in  Region  1  and  II,  and 
7 Pp  =  7 Xx  +  7 yy  =  (a  +  j/3)  p  is  the  propagation  constant  in  xy- plane.  Due  to  the  double- valued 
nature  of  the  square  root  in  the  expression  for  jzi,  mathematically  there  exist  two  types  of  solutions: 
one  has  a  positive  real  part  of  y2l  resulting  in  a  zero  field  at  infinity,  and  the  other  has  a  negative 
real  part  of  y2i  resulting  in  an  infinitely  large  field  at  infinity.  Physically,  only  the  first  group  can 
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Figure  6:  The  Poynting  vector  of  a  TE  complex  mode.  Dotted  lines  are  the  boundary  of  a  conjec¬ 
tured  box.  Solid  line  is  the  interface  of  the  dielectric  slab  and  the  air.  The  medium  and  structure 
parameters  are:  eri  =  1,  /iri  =  1,  er2  =  —2.5,  fir2  =  —2.5,  d  =  0.01  m,  k0  =  100  rad/m.  Solved 
propagation  constant:  7y//co  =  2.09195  +  jO. 86938. 


exist.  However,  this  by  no  means  implies  that  the  second  group  has  little  value.  In  fact,  the  second 
group  is  important  in  both  mathematical  manipulations  and  real  applications  as  leaky  waves  to  be 
discussed  in  detail  in  next  section  [38]. 


As  described  in  detail  in  [27]  and  listed  in  Table  1,  a  conventional  lossless  DPS  has  only  two 
types  of  proper  eigenmodes:  a  finite  number  of  discrete  ordinary  real  surface  modes,  and  continuous 
radiation  modes,  and  any  held  solutions  can  be  expressed  as  the  superposition  of  these  two  kinds  of 
modes.  With  DNG  media,  the  mode  set  becomes  much  more  complicated.  In  addition  to  discrete 
ordinary  surface  modes  and  continuous  radiation  modes,  two  discrete  surface  modes  (evanescent 
real  modes  and  complex  modes),  are  also  part  of  its  complete  mode  spectrum.  Furthermore,  unlike 
in  DPS  media  that  have  only  a  finite  number  of  discrete  surface  wave  modes,  the  DNG  media  have 
an  infinite  number  of  discrete  complex  modes  that  further  complicate  mode  solutions  [27, 39] . 
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Figure  7:  Dispersion  diagrams  for  all  modes.  Solid  lines  are  for  normalized  (3  of  the  proper  modes. 
Dashed  lines  are  for  normalized  a  of  the  proper  modes.  Dotted  lines  are  for  normalized  (3  of  the 
improper  modes.  The  medium  parameters  are:  er\  =  1,  /irl  =  1,  er2  =  —2.5,  fir 2  =  —2.5. 
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Table  1:  Complete  mode  spectrum  for  dielectric  slab  (DPS  and  DNG). 


Spectrum 

DPS 

DNG 

Discrete  evanescent  real  surface  mode 
a  =  0,  (/ 3/k0 )2  >  er2Mr2 

Improper 

Proper 

Discrete  ordinary  real  surface  mode 
a  =  0,  eri/iri  <  (/ 3/k0 )2  <  er2fir2 

Proper/Improper 

Proper/Improper 

Continuous  radiation  mode 
a  =  0,  — cx)  <  (/ 3/k0 )2  <  erinri 

Proper 

Proper 

Discrete  complex  surface  mode 
a  7^  0,  (3  7^  0 

Improper 

Proper 

3.4  Dispersion  Diagrams  and  Power  Flow 

Figure  7  shows  dispersion  diagrams  for  both  TE  and  TM  modes  for  a  DNG  slab,  including  evanes¬ 
cent,  ordinary,  and  complex  surface  modes.  Also  included  are  real  improper  modes  drawn  as  dotted 
lines.  When  the  electric  size  of  the  dielectric  slab  is  much  smaller  than  the  cntoff  electric  size  of 
the  first  real  mode,  all  complex  modes  exist  with  very  high  normalized  a  and  /3,  or  in  terms  of 
wave  properties,  very  high  attenuation  and  very  slow  phase  velocity.  As  the  k0d  increases,  f3/k0 
tends  to  decrease  rapidly  within  a  very  narrow  range  followed  by  a  steady  increase  until  its  cntoff 
point.  Notice  that  it  is  not  monotonic  and  the  value  of  f3/ko  can  be  less  than  unity  (so  called  fast 
wave),  a  notable  difference  compared  with  real  surface  modes  in  a  DPS  medium.  The  curve  of 
a/ko,  however,  monotonically  decreases  very  fast  as  k$d  increases.  At  the  cutoff  point,  a  reaches 
zero  and  (3  becomes  the  starting  point  of  a  real  mode.  The  real  surface  mode  bifurcates  into  two 
branches  from  the  cntoff  point.  One  branch  has  an  increasing  /3//c0  as  k$d  becomes  bigger  while  the 
other  has  a  decreasing  (3/k o.  The  branch  with  an  decreasing  f3/ko  will  reach  unity  shortly.  Further 
increasing  electric  size  makes  (3/ko  of  the  second  branch  begin  to  rise.  However,  it  is  no  longer 
a  proper  mode.  Note  that  the  figures  are  plotted  as  functions  of  normalized  electrical  dimension 
of  the  slab  and  we  keep  using  k0d  instead  of  frequency.  Thus  the  above  argument  can  be  used 
for  inherently  dispersive  DNG  metamaterials,  as  long  as  one  interprets  them  as  fixed  operating 
frequency  while  the  layer  thickness  is  changing. 

As  a  comparison,  the  dispersion  diagram  of  a  DPS  slab  is  plotted  in  Figure  8.  Note  that  there 
is  no  cntoff  for  TM0  mode  and  all  dispersion  curves  are  monotonic  increasing  as  the  electric  size 
increases. 

From  the  above  results,  even  for  ordinary  real  surface  modes,  DNG  media  impose  a  lot  of  dif¬ 
ficulties  by  introducing  a  bifurcation  near  its  cntoff  frequency.  A  typical  dispersion  curve  of  TE3 
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1.6 


Figure  8:  Dispersion  diagram  of  TE  modes  (dashed  lines)  and  TM  modes  (solid  lines)  for  double 
positive  material.  The  medium  parameters  are:  er\  =  1,  /ir i  =  1,  er2  =  2.5,  /iro  =  1- 


mode  is  shown  in  Fig.  9(a),  along  with  its  power  flow  in  Fig.  9(b).  From  Fig.  9(a),  it  is  seen  that 
only  a  complex  mode  exists  (branch  ‘A’)  when  kod  is  lower  than  its  cutoff  threshold.  As  shown  in 
Section  3.2,  the  zero  power  flow  in  Fig.  9(b)  shows  that  the  complex  surface  mode  does  not  carry 
away  power  in  p-direction  as  proved  in  [27].  With  the  electric  size  continuing  to  increase,  two  real 
surface  modes  show  up  from  the  cutoff  point.  The  top  branch  (branch  ‘B’)  has  a  negative  power 
flow  (power  flow  direction  is  opposite  to  the  wave  propagation  direction)  and  shows  backward  prop¬ 
erties.  When  a  waveguide  operates  at  this  mode,  its  fields  largely  are  conhned  inside  the  dielectric 
layer.  The  bottom  branch  (branch  ‘C’)  has  a  positive  power  flow  and  its  fields  extend  far  away  in 
the  air  region.  When  a  waveguide  operates  at  this  mode,  it  does  not  exhibit  backward  properties. 
Further  increasing  the  electric  size  causes  the  fields  in  the  air  region  to  decay  even  more  slowly  and 
eventually  reach  infinity.  At  that  point,  the  radiation  boundary  conditions  are  violated  and  the 
mode  becomes  improper. 

As  a  comparison,  the  TE3  mode  of  a  conventional  double  positive  medium  has  a  much  simpler 
dispersion  curve,  i.e.,  showing  up  at  the  cutoff  point  and  monotonically  approaching  the  limit  of 
\Z^r2^r2-  Readers  who  are  interested  in  detailed  derivations  and  comparisons  between  DPS  medium 
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Figure  9:  Dispersion  diagram  and  the  power  flow  in  p-direction  for  TE3  mode.  ‘A’  is  for  complex 
surface  mode;  ‘B’  is  for  top  branch  of  the  real  surface  mode;  ‘C’  is  for  bottom  branch  of  the  real 
surface  mode.  The  medium  parameters  are:  er\  =  1,  jir\  =  1,  er-2  =  —2.5,  pr 2  =  —2.5. 


and  DNG  medium  are  directed  to  [27,36]. 

3.5  The  Leaky  Modes 

Tamir  and  Kou  studied  leaky  waves  (complex  waves)  in  an  asymmetric  slab  [40]  as  shown  in  Fig¬ 
ure  10,  which  the  top  and  bottom  half-spaces  are  different.  They  claim  that  there  are  eight  different 
leaky  waves  guided  by  an  asymmetric  layer  configuration,  but  only  four  of  the  eight  can  be  excited 
in  conventional  DPS  materials.  Those  eight  leaky  waves  can  be  divided  as  four  forward  and  four 
backward  leaky  waves,  or  four  co-  and  four  contra-leaky  waves  as  shown  in  Fig.  11.  By  co-leaky 
waves  it  is  meant  that  the  fields  in  both  regions  (top  and  bottom  half-spaces)  increase  or  decrease 
exponentially.  The  contra-leaky  waves  mean  that  the  field  in  one  region  increases  exponentially, 
while  the  field  in  another  region  decreases  exponentially.  The  distance  between  the  arrows  is  in¬ 
versely  proportional  to  the  field  magnitude.  It  is  found  that  for  a  lossless  DPS  grounded  slab, 
there  are  only  two  improper  complex  modes  (leaky  waves).  If  the  slab  is  lossy,  the  conventional 
surface  wave  modes  become  complex  modes,  which  are  called  leaky  waves  by  Tamir  and  Kou.  A 
plasma  slab  (negative  e)  will  add  another  type  of  proper  complex  modes  (infinite  number  of  discrete 
complex  modes)  even  the  plasma  is  lossless.  For  a  grounded  DNG  slab  only  two  proper  complex 
modes  exist  for  a  lossless  slab.  If  the  slab  is  lossy,  the  proper  surface  wave  modes  may  become 
proper  complex  modes  or  improper  complex  modes. 
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Figure  10:  Diagram  for  an  asymmetric  slab  (top  and  bottom  half-spaces  are  different). 


Figure  12  shows  the  leaky  waves  for  a  symmetric  or  grounded  slab  with  double  positive  (DPS) 
or  double  negative  (DNG)  materials.  Tx  or  Bx  for  poles  (modes)  on  top  or  bottom  Riemann  sheet, 
x  stands  for  which  quadrant  (1st  to  4t/l).  Forward  and  backward:  forward  wave  or  backward  wave. 
It  is  found  that  for  a  lossless  DPS  grounded  slab,  there  are  only  two  improper  complex  modes  (leaky 
waves)  on  B4  and  B3.  If  the  slab  is  lossy,  the  conventional  surface  wave  modes  become  complex 
modes  on  T4,  which  are  called  as  leaky  waves  by  Tamir  and  Kou.  A  plasma  slab  (negative  e)  will 
add  another  type  of  proper  complex  modes.  For  a  DNG  slab,  there  are  proper  complex  modes  on 
T3  and  T4  for  both  lossless  and  lossy  slab. 

For  a  grounded  slab  with  DNG  materials,  only  two  proper  complex  modes  exist  for  a  lossless 
slab.  If  the  slab  is  lossy,  the  proper  surface  wave  modes  may  become  proper  complex  modes  or 
improper  complex  modes.  For  a  complete  mode  spectrum,  please  refer  to  Table  2. 

3.6  Mode  Loci  vs.  Frequency /Thickness 

Knowing  mode  distributions  and  behaviors  in  the  complex  7P-plane  is  important  in  solving  dielec¬ 
tric  waveguide  problems,  excitation  problems,  and  microwave  strip  problems.  It  provides  an  initial 
guess  when  using  Newton’s  method  to  find  solutions.  In  performing  Sommerfeld  integrations,  we 
resort  to  mode  loci  to  conceive  a  proper  integral  path  which  is  an  essential  step  in  solving  an  exci¬ 
tation  problem.  In  this  and  next  sub-sections,  we  will  first  examine  how  complex  modes  change  as 
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Forward  leaky  waves 


Directions  are  for  constant  phase  (wave  propagation).  The  source  is  at 
the  left  (z=0)  and  solutions  are  forz>0. 

_  Backward  leaky  waves 


Co-leaky  waves 


Contra-leaky  waves 


Leaky  waves  with  dashed  lines  cannot  be  excited  by  realistic  sources. 


Figure  11:  Leaky  waves  in  an  asymmetric  slab  as  shown  in  Figure  10.  The  distance  between  the 
arrows  is  inversely  proportional  to  the  held  magnitude. 


Table  2:  Complete  Spectrum  of  Dielectric  Slab 


Modes 

Finite  number, 
discrete  evanes¬ 
cent  surface 

mode 

Finite  num¬ 
ber,  discrete 
normal  sur¬ 
face  mode 

Infinite  num¬ 
ber,  discrete 
complex  mode 

Continuous 

mode 

DPS  (TE,  TM) 

Improper 

Proper/impr. 

Improper 

Proper 

Negative  e 

TE 

Improper 

None 

Improper 

Proper 

TM 

Proper 

None 

Proper 

Proper 

Negative  fi 

TE 

Proper 

None 

Proper 

Proper 

TM 

Improper 

None 

Improper 

Proper 

DNG  (TE,  TM) 

Proper 

Proper/impr. 

Proper 

Proper 

the  electric  size  k0d  varies.  Later  we  will  study  mode  loci  in  lossy  medium. 
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Directions  are  for  constant  phase  (wave  propagation).  The  source  is  at 
the  left  (z=0)  and  solutions  are  forz>0. 


B4 

T4 

T3 

B3 

Forward 

Forward 

Backward 

Backward 

DPS 

Complex 

mode 

Complex  mode 
for  lossy  slab 

Complex  mode 
for  negative  s 

Complex 

mode 

DNG 

Complex  mode 
(lossy  or  lossless) 

Complex  mode 
(lossy  or  lossless) 

For  a  given  source,  all  modes  carry  the  power  away  from  the  source. 

Figure  12:  Leaky  waves  in  a  symmetric  or  grounded  slab  for  both  double  positive  (DPS)  and  double 

negative  (DNG)  materials.  Tx  or  Bx  for  poles  (modes)  on  top  or  bottom  Riemann  sheet,  x  stands 
for  which  quadrant  (1st  to  4th).  Forward  and  backward:  forward  wave  or  backward  wave.  For  a 
DPS  slab,  there  are  improper  complex  modes  on  B4  and  B3.  If  the  slab  is  lossy,  the  proper  complex 
modes  may  exist  on  T4.  For  a  DNG  slab,  there  are  proper  complex  modes  on  T3  and  T4  for  both 
lossless  and  lossy  slab. 

For  the  most  general  case,  the  propagation  constants  are  assumed  to  be  complex  as  shown  in 
(47).  The  equations  =  —klerifiri  —  7^  (i  —  1,2  for  Region  1  or  11)  are  functions  of  7^,  therefore 
given  7P  =  a  +  j/3  a  solution,  — 7P  is  also  a  solution.  After  some  algebra  operations,  we  have  [39] 

07 1  -  Pi  1  +  2jazl(3zl  =  —k\  -  a2  +  (32  -  2 ja/3  (48) 

+  2ja22/3z2  =  ~kl  -  a2  +  (52  -  2 jaf3  (49) 

where  k'f  =  hoerj/ir?;  (i  =  1,  2).  To  satisfy  the  equations  for  the  imaginary  part,  one  must  have 

OizlPzl  =  Oiz2fiz2  =  ~Oip  (50) 
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Figure  13:  TM2  pole  loci  as  the  normalized  thickness  kod  increases.  Medium  parameters  are: 
eri  =  1,  Hr  1  =  1,  er2  =  —2.5,  p,r2  =  —1.6,  fcod  increases  from  1  to  2. 

Table  3:  Possible  sign  choices  of  complex  modes  for  lossless  medium. 


7 

7* 

-7* 

-7 

(a,  (3) 

(+>  +) 

(+> ) 

(->  +) 

ftzl) 

(+> ) 

(+>  +) 

(+>  +) 

(+> 

(aZ2,  f3z 2) 

(+,  -)  or 
(~>  +) 

(+,  +)  or 
(-  -) 

(+,  +)  or 
("  -) 

(+,  -)  or 
(~>  +) 

Simultaneously  changing  the  sign  of  f3  and  /3Z 1  does  not  violate  the  above  transcendental 
eigenequations  [39].  Therefore,  7*  is  also  a  solution.  If  7*  is  a  solution,  from  the  previous  ar¬ 
gument  we  can  conclude  that  —7*  is  a  solution  as  well.  Table  3  lists  the  possible  sign  choices 
for  a  lossless  medium.  The  first  two  columns  are  solutions  for  outward  waves  and  the  last  two 
columns  are  solutions  for  inward  waves.  Those  sign  choices  are  very  helpful  to  understand  the  be¬ 
haviors  of  complex  waves.  For  example,  for  the  outward  waves  shown  in  the  first  two  columns,  /3Z\ 
has  opposite  sign  of  f3.  The  wave  propagates  in  a  combining  +p/—z  direction  or  a  —p/+z  direction. 

When  electric  size  is  small,  only  complex  modes  exist.  According  to  the  above  conclusion, 
complex  modes  always  form  a  quartet  for  lossless  media.  Fig.  13  shows  a  particular  mode  TM2. 
With  electric  size  increasing,  complex  poles  keep  moving  toward  the  real  axis  and  become  the 
starting  point  of  the  bifurcated  real  modes.  One  of  the  two  branches  moves  inward,  touches  the 
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Figure  14:  TM2  real  pole  loci  as  loss  tangent  S  decreases  from  0.4  to  0  with  step  size  0.08.  Medium 
parameters  are:  eri  =  1,  nr i  =  1,  er2  =  —2.5(1  +  jS),  /ir2  =  —1.6,  k0d  =  1.6. 

branch  point  at  k$d  =  7t/a/3  =  1.8138  for  the  parameters  given  in  Fig.  13,  and  eventually  becomes 
improper  mode.  Meanwhile  the  other  branch  moves  outward  and  approaches  to  ^/er2/xr2^o-  This 
behavior  is  also  well-expected  from  the  dispersion  curve  in  Fig.  9(a).  The  same  comments  are  also 
applied  to  all  other  modes  shown  in  Fig.  7  except  TEi  and  TM0,  because  they  are  not  complex 
modes  for  any  electric  size. 

3.7  Mode  Loci  vs.  Loss 

In  solving  for  the  Green’s  function  of  layered  media  or  other  excitation  problems,  electromagnetic 
fields  are  often  represented  as  Sommerfeld  integrals  [31,32],  A  detailed  discussion  about  Sommerfeld 
integral  is  presented  in  next  chapter.  When  performing  integral  evaluations,  one  needs  to  specify 
the  integral  path.  There  is  no  problem  for  a  lossy  medium  because  all  the  surface  wave  poles  are  off 
the  real  axis,  and  the  inverse  Fourier  transform  along  the  real  axis  converges  due  to  the  loss.  For  a 
lossless  medium,  it  is  not  straightforward  to  conceive  a  proper  integral  path  because  the  integrand 
blows  up  near  poles  in  the  real  axis.  [31]  has  a  detailed  analysis  of  Sommerfeld  integral  path  (SIP) 
for  conventional  lossless  DPS  material.  In  short,  we  assume  a  lossy  medium  first,  then  monitor 
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Figure  15:  TM2  complex  pole  loci  as  loss  tangent  S  increases  from  0  to  0.4  with  step  size  0.08. 
Medium  parameters  are:  erl  =  1,  /url  =  1,  er2  =  —2.5(1  +  jS),  /ir2  =  —1.6,  k0d  =  1. 

how  poles  change  as  loss  becomes  smaller  and  smaller,  and  finally  construct  an  integral  path  that 
does  not  encounter  any  poles,  ffowever,  when  dealing  with  DNG  material,  we  cannot  use  exactly 
the  same  SIP  as  we  do  for  conventional  DPS  material.  Attention  has  to  be  paid  to  how  different 
kinds  of  poles  change  as  the  DNG  medium  loss  approaches  zero. 

Fig.  14  shows  TM2  pole  loci  when  loss  becomes  smaller  and  smaller.  The  blue  arrows  show  the 
direction  of  decreasing  loss.  There  are  two  kinds  of  surface  wave  modes.  The  first  one  has  most  of 
its  power  flow  in  the  air  while  the  other  has  most  of  its  power  flow  in  the  dielectric  layer.  When  the 
medium  is  lossy,  the  first  kind  of  surface  wave  poles  are  in  the  second  and  fourth  quadrants,  and 
the  second  kind  of  surface  wave  poles  are  in  the  first  and  third  quadrants.  As  the  loss  approaches 
zero,  both  merge  to  real  axis. 

Fig.  15  shows  complex  pole  loci  versus  the  medium  loss.  The  blue  arrows  show  the  direction 
of  increasing  loss.  There  are  several  observations:  First,  unlike  the  lossless  medium  where  7,  —7, 
7*,  and  —7*  are  all  solutions,  only  7  and  —7  are  solutions.  Second,  as  the  loss  becomes  larger, 
one  pair  of  complex  poles  become  faster  and  attenuate  less  (decreasing  f3  and  a),  and  the  others 
become  slower  and  attenuate  more  (increasing  (3  and  a).  This  information  may  used  in  designing 
leaky  wave  antenna  with  a  large  aperture  with  relative  uniform  held  distributions. 


4  Sommerfeld  Integral  Path  and  Numerical  Considerations 


(a)  Horizontal  magnetic  dipole  (b)  Horizontal  electric  dipole 

Figure  16:  Horizontal  magnetic  dipole  on  PEC  covered  with  DNG  media  and  horizontal  electric 
dipole  at  the  interface  of  grounded  DNG  media  and  air. 

So  far  we  have  studied  eigenproblems  of  grounded  dielectric  slabs  of  DNG  media.  When  a  finite 
source  such  as  shown  in  Fig.  16  is  considered,  the  spectral  domain  approach  is  widely  used.  At 
the  very  beginning  of  this  section,  we  will  introduce  spectral  domain  approach.  The  HMD  case 
(configured  as  Fig.  16(a))  is  solved  as  an  example  of  this  approach.  It  is  found  that  the  solutions 
are  expressed  via  a  Green’s  function.  In  a  layered  medium,  the  Green’s  function  is  written  as  a 
Sommer fcld-type  integral  [31],  whose  integrand  has  poles  on  the  fcp-plane  that  correspond  to  the 
eigenmode  solutions  of  the  same  structure  with  source  removed.  In  the  rest  of  this  section,  a  de¬ 
tailed  discussion  about  Sommerfeld  Integral  Path  is  presented. 


4.1  Spectral  Domain  Approach 

4.1.1  Hertz  Potentials 

The  Hertz  potentials  are  defined  as 

Eh  =  -ju/iW  x  nh  (51) 

He  =  jueV  x  ne  (52) 

In  source  free  and  homogeneous  isotropic  region,  both  n,  and  n/t  satisfy  the  homogeneous  vector 
Hclmholz  equation 

(V2  +  kf  )  n*’h  =  0  (53) 
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In  a  layered  structure  whose  material  constants  are  only  functions  of  z,  the  electromagnetic 
fields  are  decoupled  to  TE2  (IP  =  0.  Eh'  =  II^z)  and  TM2  (IIe  =  ne2z,  II/l  =  0)  modes  [31].  Thus 
Equation  (53)  is  reduced  to  homogeneous  scalar  Helmholz  equation  whose  cylindrical  form  is 


~ld_  d_  1  d 2 

pdp  ^  dp  p2  dcf)2 


dz 2 


j-^ez,hz 


o 


(54) 


4.1.2  Fourier  Transforms 

Given  that  the  Hertz  potentials  are  solved,  the  electric  and  magnetic  fields  are  found  by 

E  =  vxvxir-  jupv  x  n/l  =  k2ue  +  vv  ■  rr  -  jupv  x  nh  (55) 

H  =  VxVxtf  +  jueV  x  IP  =  k2Uh  +  V  V  •  Ft'1  +  jueV  x  IP  (56) 


Taking  the  advantage  of  the  symmetry  of  the  problem,  the  Fourier  transforms  adopted  by  [43] 
are  used.  They  are 


n  62(p,<M 

uhz(P,<p,z) 


—2jn  sin 
— 2j7T  cos 


Fe( A2,  z)Ji(k0Xp)dX 


P/i(A2,  z)Ji(k0Xp)d,X 


(57) 

(58) 


With  such  Fourier  transforms,  (54)  is  converted  to  1-D  differential  equation 


d2  2 
d^~lzi 


Fe,h{X,z) 


0 


(59) 


where  yA  =  kyX2  —  k2  =  k^X2  —  eripri).  Notice  that  A  =  kp/ko  is  normalized  transverse  wavenumber. 
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4.1.3  TE  Fields  in  Cylindrical  Coordinate 

The  fields  of  TE  modes  in  spectrum  domain  are 


„  .  10n'“  „  ■  j.  2  sMk o\p)^ 

Ep  =  —juju — =  2nuj  asm  cl)  /  bh  (A  ,z) - a  X 

P  <90  Jo  P 

V  ■  9Uhz  O  Jrnram 

b(j)  =  j ujp— —  =  znuj/icoscp  /  bh{ A  ,  z)  koXJ0(k0Xp) - ctA 

op  Jo  Y  P 

Ez  =  0 

„  n.  ,  r  J,  w„  ,  N 

Hp  =  =  2j7r cos cj)  /  Fft(A  ,2)  k0XJ0{k0Xp) - dA 

opoz  Jo  Y  P  J 

».  i^2nhz  n.  .  ,  r  P/,2  ^ Ji(/c0Ap) 

hp= -nnr  =  2J7rsmJ)  Fk(x  ’z) - 

p  d$dz  Jo  p 


/»oo  r  _. 

+  fc2)nh2  =  —  2  j  tv  cos  <fi  /  E"(A2,^)  +  k2Fh(X2,z )  J1(k0Xp)dX 


—  —  2jnkl  coscj)  /  X2Fh(X2,z)Ji(k0Xp)dX 
Jo 

where  i^( A2,z)  and  F[{(X2,z)  are  the  first  and  second  order  derivatives  with  respect  to  z. 


(60b) 


(60d) 


In  deriving  the  above  formula,  the  following  identity  is  used 

^dp^  =  p  [^oApdo(^oAp)  -  Ji(k0Xp)] 

When  p  =  0,  we  have  limp^o  =  JlF’x,J'1  =  the  above  formulas  are  simplified  to 

P°°  ~  k  X 

Ep  =  2mupsin0  /  Fh(X2 ,  z)-^—dX 
Jo  2 

/oo  ^ 

Fh{  A2,  z) - dX 


(62b) 


Ez  =  0 


Hp  =  —  2j7r  cos  cp  J  F'h{X2 1  z)-^-dX 

H,p  =  2j7r  sin  cj)  /  F'h{X2 ,z)-^-dX 

Jo  b, 

Hz  =  0 


(62d) 
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4.1.4  TM  Fields  in  Cylindrical  Coordinate 

The  fields  of  TM  modes  in  spectrum  domain  are 
g2]jez  _ 

E»  =  -drfI  =  -2jns'm4,l  F'(A  'z) 


k0XJ0(k0Xp )  - 


Ji(k0\p) 

P 


dX 


z,  ld2^ez  n.  J%fn2  ,Ji(koXp)^ 

E<P  =  ~^tk~  =  -2jvrcos0  /  Fe(A  ,  z) - dA 

p  d(j>dz  J0  p 


(63a) 

(63b) 


Q2 

=  (^2  +  fc2)neZ  =  -2jvrsin0  J 


F'J{ A2,  z)  +  k2Fe(X2,  z )  J1(k0Xp)dX 


—  —  2jirkl  sin  (j>  /  X2 Fe(X2 ,  z)Ji(k0Xp)dX 


u  ■  ---  „  Tsn  2 

Hp  =  jue — — —  =  Z7ru;e  cos <p  /  Fe(A  ,  z) - dA 

"  ^  70  p 

Hfj,  =  -ju:e'EJr-  =  — 2muesin0  [  Fe(X2,z )  k0XJ0(k0Xp)  -  ’h^Ml 

Jo  [  P 


ldWz 
p  d<J 
dWz 
dp 


dX 


Hz  =  0 

where  FP(X2,z)  and  Fj(X2,z )  are  the  first  and  second  order  derivatives  with  respect  to  z. 
When  p  =  0,  we  have 

/°0  _  U  \ 

K(^2,z)~Yd^ 


(63c) 

(63d) 

(63e) 

(63f) 


Es  =  —  2j7r  cos  ( 


FA  =  0 


F'e(\2 ,  z)PJd\ 


kn  X 


FIP  =  27rcae cos <p  /  Fe(X2,z)——dX 
Jo  2 

p°o  _  L  \ 

=  —  2mue  sin  0  /  Fe(A2,  z)-^—  dA 

Jo  2 

dC  =  o 


(64a) 

(64b) 

(64c) 

(64d) 

(64e) 

(64f) 


4.1.5  Rewriting  as  Infinite  Integrals 

Equation  (60)  and  (63)  are  integrals  from  0  to  oo.  It  is  cumbersome  to  do  steepest  descent  trans¬ 
formation  for  semi-infinite  integrals.  Thus  we  want  to  rewrite  them  to  integrals  from  —  oo  to  oo  by 
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using  the  formula  of 


4W  =  i[ff<1)W  +  ff<2)W] 

HW(ze~*j)  =  -evniH^\z) 


The  general  form  of  semi-infinite  integrals  involving  Jq  are 


I0=  /  4>(A2)AJ0(A)dA 

Jo 

i  roo 


*  JO 

2  TOO 


$(A2)A[^1)(A)  +  i/(S2)(A)]rfA 


2 

1 

2 


'o 

1  f00 


^>(A2)A[— i?o  (— A)  +  #<3  (A)]dA 
<f>(A2)Ai/Q2)(A)(iA 


The  general  form  of  scmi-inhnite  integrals  involving  J\  are 


/,  =  /  $(A2)J!(A)hA 

Jo 


<f>(A2)  (A)  +  if{2)  (A)]  dA 


^  Jo 

1  r°° 

2  Jo 

1  l00 


$(A2)[hfl2)(-A)  +  i/j2,(A)]dA 
$(A2)tf{2)(A)dA 


(65) 

(66) 


(67) 


(68) 


4.1.6  Total  Fields  in  Cylindrical  Coordinate 

From  the  previous  two  sections,  the  total  fields  are  written  below.  It  will  be  seen  later  that  a 
numerical  difficulty  is  avoided  by  doing  this  combination. 
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Ep  =  it  sin  < 


ufxFh(X^z)+jF^,z) 


7T  sin  0  I  —jkoXF'e(X2,z)HQ2\ko\p)dX 


H[2\k0Xp) 


P 


Ea,  =  'll  COS  ( 


upFh{X\z)+]F'e{X\z) 


H[2\k0Xp ) 


dX 


P 


/•UVJ 

upkoXFh( A2,  z)H^2\koXp)dX 

-OO 

/oo 

— X2Fe(X2,  z)H[2\k0Xp)dX 


Hp  =  7T  cos  ( 


a;eFe(A0z)+jF'(A2,z) 


+  7rcos0  /  —jkoXF{l(X2,z)HQ2\koXp)dX 


1  #|2)(Mp) 


p 


dA 


*/  — OO 

r»oo 


=  7T  sin  (j)  J 

+  7r  sin  0 

Hz  =  jnkl  cos  0  J 
If  p  =  0,  the  total  fields  are 

Ep  =  2tt  sin  0 


ueFe(X\z)+jF'h{X\z) 


H[2\k0Xp ) 


dX 


P 


—  OO 
roo 


-Luek0XFe(X2,  z)H[2\k0Xp)dX 
-X2Fh(X2,z)H[2\k0Xp)dX 


upFh{X2,z)-3F'(X2,z) 

E#  =  27TCOS0  /  LvpFh(X2,  z)  -  jF'e(X2,  z ) 

do  L 

£0  =  0 


knX 


dX 


^d\ 

2 


Hp  =  27 r  cos  0 

id <f,  =  27T  sin  0 
id.  =  0 


a)ei00A2,z)-jF'(A2,z) 


0nA 


dA 


ueFe(X2,z)-jFi(X2,z) 


'XdX 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 
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4.1.7  Solving  Fe  and  Fh 

Now  let  us  solve  the  equation  (59)  for  Fei  and  Fhi  (i  =  1,  2).  The  air  region  (region  1)  is  unbounded, 
so  the  general  solution  of  Hertz  potential  is 

Fe  i  =  Ae-^l{z~d)  (81) 

Fhl  =  De~^l{z~d)  (82) 

The  Hertz  potential  in  the  dielectric  layer  (region  2)  is 

Fe2  =  B  sinh  [7*2 (d  -  z)\  +  C  cosh  [7 z2(d  -  z)}  (83) 

Fh 2  =  E  sinh  [7 z2(d  —  z)\  +  F  cosh  [7*2(0?  -  z)]  (84) 


At  medium-air  interface,  applying  boundary  conditions  for  TE  and  TM  modes,  we  get 

A* rlFhl  =  hr  2^2 

dFhi  _  Fh2 

dz  dz 

0  r  1  1  0  r2  i‘  t  2 

d_FA  =  F\2 
dz  dz 

Now  we  have  four  equations  with  six  unknowns 

er\A  =  er2C 

IziA  =  7 z2B 

Hr  1 D  =  HrlF 

lz\D  =  7*2  A 


(85) 

(86) 

(87) 

(88) 


(89) 

(90) 

(91) 

(92) 


We  need  two  more  equations  to  solve  the  problem.  They  are  found  by  invoking  the  boundary 
condition  fxEo  =  —  at  z  —  0  (the  tangential  component  of  the  electric  field  has  a  discontinuity 
if  there  is  a  finite  magnetic  surface  current).  But  the  TE  and  TM  modes  are  no  longer  decoupled 
anymore  at  z  —  0  and  total  electric  holds  must  be  used 

E2x  I z=0  =  E2p  cos  (f)  -  E2<t,  sin  0  =  0  (93) 

E2y\z=0  =  E2p  sin  0  +  E2^COS0  =  Pm5(x)5(y)  (94) 

where  Pm  is  the  magnetic  dipole  moment  (magnetic  current  times  length). 
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From  (93),  we  have 


The  left  hand  of  (94)  is 


+  jF'e2  —  0| 


2=0 


(95) 


E2p  sin (j)  +  E^  cos  (j)  —  2nujp2  /  k0XFh(X2 ,  z)J0(k0Xp)dX 

Jo 


(96) 


In  the  right  hand,  using  the  identity  in  [31],  p.  575,  we  have 


p  p  /*oo  p  7,2  /*oo 

PmS(x)S(y)  =  ~^5{p)  =  /  XJ0(Xp)dX  =  -Jg*  /  XJ0(k0Xp)dX 

2vrp  27t  y0  27t  y0 


(97) 


In  equating  (96)  and  (97),  one  get 


-2jirk0F'e2  =  2irup2k0Fh2 


P  k2 
1  m1^  o 

2vr 


or 


2  =  0 


Pmko 

4:7T2Up2 


Now  we  have  six  equations  and  six  unknowns.  Upon  solving  the  equations,  we  get 


A 

B 

C 

D 

E 

F 


j  koPm^r2 
Ati2Dtm 
j  k()Pm^r2'y zl 
4tt2U)TM722 
j  koPm^rl 
47 T2DtM 

Pmko'Y  z2  pr2 

4:7T2Up2DTE 

Pmko'Y  zlpr2 

4:7T2Up2DTE 

Pmko'Y  z2prl 

Att2ujp2Dte 


(98) 


(99) 
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And  for  Fe  h,  we  have 


Fe  1 
Fe2 
Fhl 
Fh2 


jPmk  per2  ^Zl(z-d) 
An2DTM 


j  Fm  k(  i 


4:TT2DtM^z2 

Pmlz2  f 

47 T27]0DTE 
P 

L  ID. 


{er2lzi  sinh[7 z2{d  -  z)\  +  eri7z2  cosh[7~2(d  -  2:)]} 


=-7zi  (z-d) 


4:7l27]0DTEldr2 


{^r2lzi  sinh[7z2(d  -  z)\  +  pri7 z2  cosh[7z2(d  -  z)]} 


where 


Dte  =  Hr2lz\  sinh(7z2d)  +  [ir\lz2  cosh(-fz2d) 

DTm  =  tr2lzi  cosh(7z2d)  +  eri7z2  sinh(7z2d) 

They  are  the  same  as  the  eigen  equations  we  derived  earlier. 

Their  derivatives  with  respect  to  z  are 

771/  _  jPmer2h'yzl  Jzl(z-d) 

el  4:7t2Dtm 

K 2  =  /Tn^0  {er2lzi  cosh[7z2(d  -  z)\  +  eri7z2  smh[jz2(d  -  z)}} 

A,kzDtm 

pi  __  PynP  z2^  z\  —7 zl(z—d ) 

hl  ~  An2Vo Dte 

P'h2  =  ~A — 9 — —  {Hr2lzi  cosh[7z2(d  -  z)\  +  /7ri7z2sinh[7 z2(d  -  z)}} 
A7l21]0DTEldr2 


(100a) 

(100b) 

(100c) 

(lOOd) 

(101) 

(102) 


(103a) 

(103b) 

(103c) 

(103d) 


4.2  Sommerfeld  Integral  Path  (SIP) 

After  the  proper  integral  path  has  been  conceived,  it  is  often  deformed  to  the  branch  cut  (BC)  to 
provide  a  spectral  representation  (combination  of  all  possible  proper  modes  including  both  discrete 
and  continuous  modes)  or  to  steepest  descent  path  (SDP)  shown  in  Fig.  17  to  approximately  the 
far  held  efficiently  [31].  According  to  Cauchy  theorem, 

I  sip  =  Inp  ~  Znj  ^2  R .es[f(kp),  kpi ]  (104) 
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where  Isip  is  the  integral  along  the  Sommerfeld  integral  path  (SIP),  I^p  is  the  integral  along  a 
new  path  (NP),  which  is  usually  chosen  to  be  a  BC  or  a  SDP.  The  last  term  is  the  contribution 
from  all  poles  enclosed  by  these  two  paths. 


□  Branch  point 

- SIP 

—i —  Branch  cut 
■  - -SDP 
O  Proper  modes 
X  Improper  modes 
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T1 


I  I  I  I  I  i  I  II 
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o .. 


ffl  - 


I  I  I  I  I  I  I 
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X 


* 


i . . 
/  .. 
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B4 


B1 


Figure  17:  Poles  and  integral  paths  in  transformed  0-plane.  Medium  parameters  are:  er2  =  —2.5  — 
jO.l,  fir2  =  —2.5,  k0d  =  2.  Observation  angle  9  =  60°  (z  =  0). 


As  one  of  the  most  commonly  used  transformations  by  mathematicians  and  physicists,  the 
conformal  mapping  of 

kp  —  k\  sin  0  (105) 

where  k\  is  the  wavenumber  in  air  region  and  (j>  =  a+jrj  is  a  complex  number,  maps  the  overlapped 
Riemann  sheet  to  an  unfolded  0-plane  [36,41].  Fig.  17  plots  poles  and  integral  paths  and  branch 
cuts  in  0-plane  in  a  much  clearer  way.  It  is  clearly  shown  that  real  surface  modes,  complex  surface 
modes  in  the  third  and  fourth  quadrants,  and  continuous  radiation  modes  from  the  branch  cuts 
form  a  complete  mode  spectrum  of  DNG  medium. 
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When  the  medium  is  lossy,  all  poles  are  not  along  the  real  axis  in  general.  So  the  Sommerfeld 
integrand  can  be  evaluated  along  the  real  axis  without  any  difficulty.  But  when  the  medium  is 
lossless,  the  poles  are  located  on  the  real  axis.  All  integrands  in  the  last  section  have  Dte  or  Dtm 
in  their  denominator.  Thus  those  integrands  go  to  infinite  at  those  poles  along  the  real  axis.  A 
common  way  to  address  this  difficulty  is  to  assume  the  medium  lossy  and  then  monitor  how  the 
pole  loci  change  as  the  loss  becomes  smaller  and  smaller.  For  the  lossy  medium,  the  SIP  is  along 
the  real  axis.  When  the  loss  is  zero,  the  poles  are  on  real  -axis  and  the  SIP  should  be  deformed 
lettle  away  from  those  poles.  From  Fig.  14,  it  is  quite  natural  that  one  chooses  the  SIP  as  shown 
in  Fig.  18(a).  As  a  comparison,  Fig.  18(b)  shows  the  SIP  for  a  DPS  medium.  From  Fig.  9,  it 
is  observed  that  the  first  kind  of  poles  corresponds  to  forward  surface  modes  (the  power  flows  in 
the  same  direction  as  the  phase  velocity),  while  the  second  kind  of  poles  corresponds  to  backward 
surface  modes.  When  the  medium  is  lossy,  all  first  kind  poles  are  located  in  the  second  and  fourth 
quadrants  and  all  second  kind  poles  are  in  the  first  and  third  quadrants.  Now,  consider  the  cases 
shown  in  Fig.  16.  The  branch  C  in  the  figure  is  the  first  kind  and  the  branch  B  is  the  second  kind. 
Because  the  power  always  leaves  from  the  source,  one  must  include  the  first  type  of  poles  with 
kp  >  0  and  the  second  type  of  poles  with  kp  <  0  as  shown  in  Fig.  18(a).  So  both  included  modes 
carry  away  power  from  the  source  even  they  have  opposite  phase  velocity. 


(a)  DNG  media  (b)  DPS  media 

Figure  18:  Poles  and  Sommerfeld  integral  paths  in  fcp-plane. 
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4.3  Residue  at  Origin 


Usually,  the  Sommerfeld  integrand  involves  Bessel/Hankel  functions  for  layered  media  [31].  As  the 
integral  variable  approaches  the  origin,  the  argument  of  the  Bessel/Hankel  function  approaches 
zero.  The  small  argument  first  order  Bessel  function  can  be  approximated  as  [42] 


Ji{z)  « 

Y^z)  « 


z  z 


2T(2)  2 

2  2 

- V  —Ji(z)  In 

7 TZ  7 r 


z 

2 


(106) 

(107) 


(2) 

Therefore,  the  first  order  Hankel  function  H\  \z)  =  J\{z)  +  j]Y\  (z)  has  a  simple  pole  at  origin 
introduced  by  its  imaginary  part  Y\  (z).  Special  care  must  be  taken  to  account  for  it  otherwise 
the  integral  can’t  be  performed.  Generally,  there  are  two  ways  to  overcome  this  difficulty.  As  an 
example,  we  will  study  a  dipole  at  the  interface  of  a  grounded  dielectric  slab  and  air. 


The  first  method  is  to  simply  combine  the  TE  fields  and  TM  fields  together  in  order  to  calculate 
the  total  field.  From  [39,43],  the  total  electrical  field  at  p-direction  reads 


EP  = 


ufiFh(X2,z)+jFi(X2,z) 


H[2\k0Xp ) 


dX 


P 


Fg(A2,  z)koXHQ2\koXp)dX  />  7r  sin  </> 


where  A  =  kp/kG  and  the  details  of  Fe  and  Fh  can  be  found  in  4.1.7. 


(108) 


There  is  no  first  order  pole  in  the  second  term.  For  the  first  term,  up>Fh(X2,  0] 

upFh(X2,z)+jF'(X2,z) 


7^  0  and 
is  zero  at 


e  at  origin  and  the  integral  will  not 


jFe'(A2,0)  7^  0.  However,  as  proved  in  [39],  their  summation 
origin.  Therefore,  we  have  a  zero  to  cancel  out  the  simple  po’ 
encounter  any  problems. 

When  a  combination  of  the  TE  and  TM  fields  is  not  allowed,  such  as  when  one  wants  to 
calculate  the  contributions  from  TE  and  TM  modes  separately,  we  have  to  resort  to  the  second 
method:  treat  TE  and  TM  case  separately.  Usually,  it  can  be  proved  to  be  a  Cauchy  principal 
value  (CPV)  problem  and  calculated  as 

pO—e  po 


I  =  lim 

e— >0 


+ 


'  0+e. 


f{z)dz 


(109) 


It  is  more  complicated  when  one  deforms  the  integral  path  to  a  BC  or  a  SDP  as  mentioned  pre¬ 
viously.  In  this  case,  one  has  to  artificially  add  a  small  segment  to  the  integral  path  to  close  the 
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Figure  19:  Integral  paths  near  the  origin. 


integral  path  since  we  cannot  apply  Cauchy’s  theorem  to  a  CPV  integral  (because  it  has  a  gap), 
and  later  subtract  it  from  the  result.  As  an  example,  let  us  take  a  look  at  the  TE  component  of 
equation  (108) 


Ep  =  7ru;/isin0  /  Fh(X2,z) 


Hi2\k0Xp) 


dX 


'  —  OO 

poo 


p 


2  ^Ji(k0Xp) 


—  2n(jjp sin (f>  /  Fh(X2,z) 

Jo  P 

Because  the  Bessel  function  Ju{z)  is  an  entire  function,  we  have 


dX 


(110) 

(111) 


h  = 


f(kn)Ju(kpP)dkp  =  lim  /  f(kp)Ju(kpp)dkp 


£— >0 


Now  if  written  in  Hankel  form,  the  integral  I o  is 

rO—e  pc 


In  =  lim 


e — >0 


—  oo  J  0+e_ 


f(kr)H^\kpf,)dkf 


(112) 


(113) 


Therefore,  the  original  integral  (110)  is  actually  a  CPV  problem  with  integral  path  shown  as  a  solid 
black  line  in  Fig.  19. 

When  the  integral  path  is  detoured  to  a  BC,  in  order  to  apply  the  Cauchy  residue  theorem,  one 
has  to  cover  the  infinitesimal  small  gap  of  the  CPV  path  by  deliberately  adding  a  small  segment 
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at  origin  as  drawn  with  a  dashed  line  in  Fig.  19.  Now,  with  that  small  segment  added,  one  has 

Isiph  =  I  sip  +  7 rjR+  (114) 

where  Isip  is  the  original  integral,  Isiph  is  the  integral  with  the  small  gap  covered,  and  R+  is  the 
residue  at  the  positive  edge  of  the  branch  cut.  Similarly  for  BC  path, 

Ibch  =  Ibc  +  (R+  —  R  )  (115) 

where  Ibc  is  the  CPV  integral  along  BC,  Ibch  is  the  integral  with  the  small  gap  covered,  and  R+ 
and  R~  is  the  residue  at  the  positive  and  negative  edges  of  the  branch  cut.  Applying  the  Cauchy 
residue  theorem  to  Isiph  and  Ibch ,  one  has 

Isiph  =  Ibch  —  27T j  Res[/(z),  k,]  (116) 


Finally,  one  arrives  at 


Isip  =  Ibc  -  ^(R+  +  R~)  -  2 irj  £  Res \f(z),  kt] 

i 

The  calculation  of  the  residue  at  the  origin  is  straightforward.  When  A 
approximated  as 


(117) 

0,  the  integrand  is 


Fh(\\z) 


2  ^H[2\k0Xp)  Fh(0,z) 


P  P 

„  Fh(0,z)  2 j 


[Ji(k0\p)  -jYi(koXp)] 


p  irkoXp 


(118) 


The  residue  of  the  integrand  at  the  origin  is  then 


R± 


2jjUM 

ukop2 


7zl=±?fcl 


(119) 


Fortunately,  it  is  only  in  rare  cases  that  one  needs  to  pay  attention  to  the  aforementioned 
details.  Most  of  the  numerical  calculations  can  be  performed  using  the  Bessel  function  form  of  the 
Sommcrfcld  integrals.  Therefore,  the  original  [— oo,  oo]  integral  range  is  folded  to  [0,  oo],  and  there 
is  no  problem  at  the  origin.  One  must  pay  attention  to  including  the  pole  at  the  origin  only  when 
the  complex  plane  analysis  is  needed. 
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Convergence  Test:  Hz  kernel,  A  :  0  — >  1 


Figure  20:  Convergence  tests  of  equation  (120)  and  (123). 


4.4  Numerical  Integration  Along  SIP 

It  is  found  that  the  convergence  of  the  Sommerfeld  integral  along  SIP  in  fcp-plane  is  very  slow.  As 
an  example,  we  test  two  integrals  that  have  exact  solutions.  They  are  Hz  and  Ez  in  free  space 


Hz 

Ez 


X2e~lzlZ  Ji(k0\p)d\ 


X2—e-'"lsJ1(k0\p)d\ 

721 


(120) 

(121) 


For  simplicity  purpose,  both  Hz  and  Ez  in  above  equations  are  normalized  with  the  coefficient: 


— 2j7r/e(jPm  cos  (f) 
4vr2?/ohri 


(122) 


In  (120),  the  integrand  is  smooth  along  the  integral  path.  However,  its  first  derivative  has  singu¬ 
larity  at  A  =  1.  For  this  reason,  the  convergence  is  very  slow.  The  integrand  of  (121)  is  even  worse 
as  its  integrand  blows  up  at  A  =  1.  But  it  is  still  integrable  at  A  =  1. 
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The  above  two  equations  are  transformed  to  A2i-plane,  and  now  A2i  =  kzi/k0  becomes  the 
integration  dummy  variable.  Also  changed  is  the  integral  path  that  is  mapped  to  a  different  shape 
in  A2i-plane 


Hz 

Ez 


AA2ie  7zl" Ji(/coAp)dA2i 

Ae_7zl2Ji(/c0Ap)(iA2i 


(123) 

(124) 


The  comparisons  of  convergences  of  (120)  and  (123)  are  plotted  in  Fig.  20.  The  integration  in  /em¬ 
plane  (or  A2i-plane)  is  far  superior  to  the  integration  in  Ay-plane  (or  A- plane).  Similar  conclusions 
apply  to  (121)  and  (124).  The  reason  is  because  (123)  and  (124)  now  have  continuous  high  order 
derivatives  and  thus  their  convergence  is  much  faster  than  (120)  and  (121). 
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5  Numerical  Examples  and  Optimization  for  Directivity 

In  this  section,  two  important  radiation  problems  —  a  horizontal  magnetic  dipole  (HMD)  on  PEC 
covered  with  DNG  media  and  a  horizontal  electric  dipole  (HED)  at  the  interface  of  grounded  DNG 
media  and  air  shown  in  Fig.  16  —  are  solved  using  Sommerfeld  integrations.  They  are  of  great 
importance  in  that  they  provide  Green’s  function  to  slot  antennas  and  microstrip  patch  antennas, 
because  Green’s  functions  are  electric  and  magnetic  fields  generated  by  electric  dipoles  on  the 
interface  or  magnetic  dipoles  on  the  metal  surface.  We  also  investigated  the  radiation  intensity  and 
directivity  of  electric/magnetic  dipoles  over  a  grounded  slab  to  find  optimal  values  of  the  material 
parameters  and  thickness  for  a  given  frequency. 

5.1  Numerical  Example  of  New  SIP  Path  for  DNG  Media 


Figure  21:  The  examination  of  the  proper  Sommerfeld  integral  path  for  DNG  medium,  k^d  =  2, 
er2  =  —2.5  —  jd,  nr2  =  —2.5.  Solid  line  is  for  5  =  0.01;  dashed  line  is  for  5  =  0  and  our  proposed 
SIP;  dash  dot  line  is  for  6  =  0  and  conventional  SIP. 

The  first  numerical  example  is  to  verify  the  validation  of  the  newly  proposed  Sommerfeld  integral 
path  for  DNG  media.  The  electric  fields  of  a  horizontal  electric  dipole  at  the  interface  of  air  and 
a  grounded  dielectric  DNG  medium  slab  are  calculated  along  the  interface  by  our  program  which 
can  handle  both  lossy  and  lossless  media.  In  Fig.  21,  the  solid  line  is  for  a  lossy  medium  whose 
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permittivity  is  er2  =  — 2.5  —  O.Olj.  As  the  loss  keeps  all  surface  mode  poles  off  the  real  axis,  i.e.,  kp- 
axis,  we  do  not  encounter  any  numerical  troubles.  More  importantly,  the  SIP  can  be  unambiguously 
defined  along  the  real  axis  and  both  ordinary  surface  wave  modes  with  the  first  kind  of  pole  (such  as 
TM2  )  and  with  both  kind  of  pole  (such  as  TE3)  are  excited,  which  can  be  confirmed  by  Fig.  7  with 
k0d  =  2.  The  dashed  line  in  Fig.  21  shows  the  total  field  (as  a  function  of  distance  from  source)  for 
the  same  medium  with  loss  removed,  using  the  new  SIP  proposed  in  Section  4.  It  is  seen  that  these 
two  lines  are  closely  lined  up  at  the  near  field  where  the  attenuation  is  hardly  appreciable.  As  we 
move  away  to  the  far  field,  the  electric  field  in  the  lossy  medium  undergoes  more  attenuation  and 
we  see  more  deviations  from  the  dashed  line,  but  they  exhibit  the  same  oscillation  behavior  with 
the  same  rhythms.  This  proves  the  correctness  of  the  new  SIP.  The  dash-dot  line  is  the  calculation 
for  the  same  lossless  medium,  except  that  the  SIP  for  DPS  medium  is  used,  as  shown  in  Fig.  18(b). 
It  deviates  from  either  the  blue  solid  line  or  dashed  line.  In  summary,  Fig.  21  justifies  our  proposed 
Sommerfcld  integral  path. 

5.2  The  Discrete  Pole  Contribution  of  DNG  Media 

The  second  example  solves  a  horizontal  magnetic  dipole  on  a  PEC  covered  by  a  DNG  slab  as 
shown  in  Fig.  16(a).  The  problem  is  also  studied  by  P.  Baecarelli  et  al.  [18].  The  permittivity  and 
permeability  of  this  material  are  modeled  as  [6,2] 

^  [  2  2 

pr(  u)  =  1 - - - ^  and  er(cn)  =  1 - §  (125) 

uA  —  (Uq  uA 

where  F  =  0.56,  ujq/2i t  =  4  GHz,  and  ujp/2h  =  10  GHz.  The  material  simultaneously  exhibits 
negative  permittivity  and  permeability  from  /  =  4  GHz  to  6.03  GHz  [18]. 

Figure  22  plots  the  total  electric  field  as  a  blue  solid  line  and  it  exactly  agrees  with  [18].  Also 
plotted  in  red  dashed  line  is  the  electric  field  from  the  discrete  mode  contribution.  This  includes 
the  evanescent  surface  mode  TMo  (first  kind),  complex  surface  modes  TM2  and  TE3,  and  all  other 
higher  order  complex  modes.  A  very  interesting  phenomenon  is  observed:  unlike  a  DPS  medium, 
which  always  has  at  least  one  surface  pole  whose  contribution  goes  to  infinity  as  the  observation 
point  moves  towards  the  source,  the  DNG  medium’s  large  number  of  discrete  poles  tends  to  have 
some  cancellation  effect  that  limits  their  contribution  to  a  finite  value  even  when  the  observation 
point  is  very  close  to  the  source.  Notice  that  in  our  calculation,  400  modes  are  used  in  order  to 
achieve  accurate  pole  contributions.  For  the  field  of  single  mode,  it  should  be  inverse  proportional 
to  the  square  root  of  p,  But  our  numerical  results  from  the  figure  show  that  the  total  pole  con¬ 
tributions  (with  400  poles)  become  finite  when  the  p  approaches  zero  (the  red  dashed  line).  The 
similar  plot  is  also  given  in  [18]  with  smaller  number  of  modes  included. 
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Figure  22:  The  comparison  of  our  results  to  reference  [18].  The  medium  parameters  are  exactly 
the  same  as  in  [18],  the  thickness  is  60  mm,  and  the  frequency  is  5.6  GHz. 


Figure  23:  The  electric  fields  due  to  the  pole  contributions  and  the  difference  between  the  integrals 
along  the  SIP  and  BC.  The  medium  parameters  are  given  by  (125),  the  frequency  is  4.8  GHz,  and 
the  thickness  is  20  mm. 
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In  Figure  23,  we  plot  the  electric  fields  due  to  the  pole  contributions  and  difference  between  the 
integrals  along  the  proposed  SIP  and  along  the  BC.  The  frequency  is  4.8  GHz  and  the  thickness 
of  the  slab  is  20  mm.  At  this  electric  size,  the  evanescent  surface  mode  TEi  (second  kind)  and 
ordinary  surface  mode  TM2  (both  kinds),  and  complex  surface  modes  are  excited  [21].  Excellent 
agreement  is  observed  from  Fig.  23  and  the  relation  shown  in  (104)  is  satisfied,  which  verifies  the 
correctness  of  our  new  proposed  Sommerfeld  integral  path  for  DNG  medium. 

5.3  3-D  Radiation  Patterns  (linear  scale)  of  HED  and  HMD  with  DNG 
Media 

The  radiation  patterns  associated  with  the  structures  in  Fig.  16  are  shown  in  Fig.  24.  The  media 
parameter  defined  in  (125)  is  used  in  the  calculations.  The  frequency  is  chosen  at  /  =  5.943  GHz 
as  it  gives  the  best  directivity  [18,39].  The  medium  parameters  at  this  frequency  are  er2  =  —1.831 
and  nr2  =  —0.0238.  The  thickness  of  the  dielectric  slab  is  d  =  60  mm  ( k^d  =  7.47  and  kd  =  1.56). 
Maximum  directivities  are  observed  at  the  ^-direction  (broadside):  26.8  dB  for  HMD  and  a  smaller 
number  17.2  dB  for  HED  due  to  some  cancellation  of  the  image  from  the  ground.  They  are 
drastically  larger  than  the  values  of  8.5  dB  and  12.0  dB,  respectively,  obtained  for  a  conventional 
DPS  slab  with  (er  =  2.55  and  fir  —  1  at  /  =  5.4  GHz).  Unlike  the  electromagnetic  fields  in 
conventional  DPS  media  where  some  of  the  energy  is  carried  away  by  surface  wave,  the  fields  in 
DNG  medium  tend  to  decay  more  slowly  due  to  the  existence  of  complex  modes.  Therefore,  the 
dipole  with  DNG  media  has  a  much  larger  equivalent  aperture  size  that  results  in  a  much  narrower 
radiation  beam.  The  microstrip  antennas  with  DNG  media  can  have  compact  size  comparing  the 
antennas  with  conventional  DPS  meida. 

5.4  Optimization  for  Directivity 

One  of  the  research  interests  is  to  find  optimized  values  for  the  constitutive  parameters  of  the  struc¬ 
ture  shown  in  Fig.  16,  for  a  specified  k^d  so  that  we  can  obtain  high  directivity  in  the  broadside 
direction.  The  HMD  case  (Fig.  16(a))  is  taken  as  an  example  in  the  following  discussion.  In  refer¬ 
ence  [18],  it  was  found  that  that  highly  directive  beams  can  be  obtained  by  employing  substrates 
with  a  large  ratio  of  |er|/|/xr|,  because  the  broadside  power  density  P{6  =  0)  is  proportional  to  this 
ratio. 

We  further  optimize  kod,  er,  and  /ir  to  get  the  maximum  directivity.  Figure  25  shows  how  the 
broadside  directivity  changes  along  er  and  /ir  for  a  given  k0d: 

Figure  25  gives  us  an  idea  that  for  a  given  k0d,  the  broadside  directivity  may  reaches  a  maxi¬ 
mum  when  the  product  of  er  and  /ir  is  equal  to  a  certain  value.  This  idea  is  confirmed  by  a  further 
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(b)  Radiation  pattern  for  Fig.  16(b) 

Figure  24:  Radiation  patterns  of  HMD  and  HED  associated  with  DNG  media  at  /  =  5.943  GHz. 


investigation  of  the  data  constructing  Fig.  25.  From  that  data  we  can  also  get  the  value  of  er/iT  for 
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Figure  26:  Broadside  directivity  as  a  function  of  er  (ernr  =  0.062,  dj A0  =  1). 

dj Ao  =  1  ( k0d  =  2n)  is  about  0.062.  Numerical  results  of  cases  with  different  k0d  values  lead  to 
the  same  conclusion:  for  a  given  k0d,  the  broadside  directivity  is  only  a  function  of  er.  Figure  26 
is  presented  as  a  validation  of  this  conclusion.  Figure  26  also  shows  that  a  larger  \er\  gives  a  larger 
broadside  directivity.  Because  the  product  of  er  and  /ir  is  kept  constant,  the  larger  er  is,  the  larger 
is  the  ratio  of  |er|/|//r|.  This  agrees  with  the  finding  in  [18]. 

To  verify  that  er/zr  is  unique  to  have  maximum  directivity  for  a  given  d/ Ao,  we  examine  for  a 
small  range  of  er/ir.  Figure  27  is  the  contour  plot  of  maximum  broadside  directivity  as  a  function 
of  er  and  the  product  of  er/ir.  From  this  figure  we  can  see  that  for  any  er,  the  maximum  directivity 
is  always  obtained  when  er/ir  is  the  same  certain  value  (some  value  between  0.06  and  0.065).  So 
we  can  conclude  that  to  have  maximum  broadside  directivity  for  a  given  d/ Ao,  the  product  of  er/ir 
should  be  a  unique  value  c(k0d )  =  erfir.  Hereby  we  investigate  the  dependence  between  c(k0d )  and 
electric  size.  For  the  electric  size  in  the  range  of  our  interest,  that  is  0.1  to  1  wavelength,  we  find 
the  optimized  c(kod)  at  10  different  electric  sizes  between  0.1  and  1  wavelength  and  fit  the  data 
with  a  simple  curve.  The  fitting  result  is  presented  in  Fig.  28  and  the  fitted  formula  is  given  as 

„  A  2.4635  0.0624 

c(  0 } "  I W  “ 
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where  c(k0d )  =  er/ir. 


Figure  27:  Contour  plot  of  the  directivity  at  the  broad  side  over  different  er  and  er/ir,  dj Ao  =  1. 

Now  for  any  given  electric  size  [dj Ao  or  k0d)  and  relative  permittivity  er  we  can  calculate  the 
optimized  relative  permeability  //r ,  so  that  we  can  get  the  maximum  broadside  directivity. 


From  the  discussion  above,  we  can  summarize  a  procedure  to  get  optimized  constitutive  para¬ 
meters  for  a  given  hod: 

1.  For  any  given  d/X0  (or  k0d),  we  set  er  to  a  proper  value.  Any  er  in  the  range  specified  by  the 
application  will  do;  it’s  not  important  in  this  step  because  we  just  want  to  find  the  product  of 
er  and  /ir  here,  which  is  er  independent.  And  then  we  do  a  one  variable  optimization  to  find 
the  product  of  er  and  /ir  so  that  the  maximum  broadside  directivity  is  obtained.  In  Matlab, 
there  are  certain  commands  which  can  help  with  the  optimization. 

2.  If  there  is  a  range  limit  for  er,  the  directivity  at  the  broad  side  is  always  obtained  at  the  largest 
absolute  value  of  er.  The  corresponding  value  of  /ir  is  obtained  by  using  /ir  =  c(k0d)/er. 

From  the  discussion  above,  we  know  that  the  maximum  broadside  directivity  can  be  expressed 
as  a  function  of  electric  size  d/Xo  and  relative  permittivity  er.  Now  we  provide  an  approximate 
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Figure  28:  c(k0d )  as  a  function  of  electric  size. 


formula  to  estimate  the  maximum  broadside  directivity  of  this  structure.  To  get  the  formula,  the 
maximum  broadside  directivities  for  different  d/ Ao  and  er  are  calculated  numerically.  We  use  a 
formula  with  simple  form  and  nn-determined  coefficients  to  fit  the  data  and  get  the  coefficients. 
Finally  we  find  that  the  maximum  broadside  directivity  D0  can  be  obtained  using  the  approximate 
formula: 


D0  = 


(127) 


Figure  29  shows  the  validation  of  our  approximate  formula.  By  using  Equation  (127),  we  can 
obtain  an  estimated  maximum  broadside  directivity  for  a  given  electric  size  and  relative  permittiv¬ 
ity.  This  formula  is  derived  under  the  condition  that  \er\  »  /ir |  and  for  a  given  er  and  /ir  that 
satisfy  Equation  (126).  The  numerical  results  in  Fig.  29  show  that  this  formula  works  in  the  range 
0.2  <  d/ A0  <  1. 
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Figure  29:  Approximate  formula  (Equation  (127))  vs  numerical  results. 


6  Conclusions,  Publications,  and  Significant  Events 

In  this  final  report,  we  first  review  the  mode  properties  of  a  grounded  DNG  medium  slab.  The  com¬ 
plexity  of  the  mode  spectrum  and  pole  behavior  under  parameter  and  frequency  variations  demands 
a  new  Sommerfeld  integral  path  that  is  different  from  the  SIP  for  a  conventional  DPS  medium, as 
shown  in  Fig.  18(b).  Two  numerical  issues  associated  with  Sommerfeld  integrals  are  studied  in 
great  detail.  A  new  Sommerfeld  integral  path  is  suggested  to  increase  numerical  convergence  by 
eliminating  singularities  of  the  integrand  or  the  derivatives  of  the  integrand.  Numerical  examples 
are  given  to  validate  the  necessity  of  a  new  Sommerfeld  integral  path.  The  electromagnetic  fields 
and  radiation  patterns  of  horizontal  electric  and  magnetic  dipoles  with  slabs  of  DNG  media  are 
plotted.  Due  to  the  slowly  decaying  fields  along  the  interface,  the  radiation  directivity  is  greatly 
enhanced. 

We  investigated  the  radiation  intensity  and  directivity  of  electric/magnetic  dipoles  over  a 
grounded  slab  to  find  optimal  values  of  material  parameters  and  thickness  for  a  given  frequency. 
We  optimize  the  directivity  over  the  thickness  and  the  constitutive  parameters.  It  is  found  that  for 
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a  slab  with  a  fixed  electric  size  (proportional  to  the  product  of  the  thickness  and  frequency),  the 
directivity  reaches  the  maximum  when  the  product  of  er  and  /ir  is  constant,  which  depends  on  the 
electric  size  only.  The  maximum  directivity  is  proportional  to  the  ratio  between  er  and  /J,r.  Finally, 
the  optimized  structure  to  get  maximum  directivity  at  broadside  is  investigated  and  fitted  formula 
for  the  maximum  directivity  is  provided. 


In  summary,  we  have 

1.  Found  complete  spectrum  in  ground  DNG  slab,  studied  surface  wave  suppression. 

2.  Derived  the  power  flow  for  all  discrete  surface  wave  modes  and  found  that  the  complex  wave 
modes  in  lossless  slab  do  not  carry  away  any  power. 

3.  Plotted  loci  of  complex  (leaky)  wave  modes  as  functions  of  thickness  and  loss. 

4.  Derived  field  expressions  for  a  HMD  and  HED  over  a  grounded  DNG  slab  and  found  that  the 
radtion  patterns  have  large  directivity. 

5.  Developed  empirical  formula  for  the  directivity  at  broadside. 

6.  Found  the  correct  Sommerfeld  integral  path  for  a  DNG  slab. 
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Publications 


W.W.  Shu  and  J.M.  Song,  ’’Sommerfeld  Integral  Path  for  Layered  Double  Negative  Metamate¬ 
rials,”  IEEE  Transactions  on  Antennas  and  Propagation,  under  review. 

Collaborations 

We  have  collaborated  with  Dr.  Shore  (AFRL  Hansom)  and  Dr.  Yaghjian  (self-employed)  on 
leaky  wave  properties  on  dielectric  slab  and  hybridization  of  analytical  and  numerical  approaches 
to  study  the  band  structures  of  metamaterials. 

Significant  Events 

1.  The  PI  attended  the  progress  review  meeting  on  April  27-28,  2010  in  Arlington,  Virginia,  and 
briefed  the  audience  on  our  latest  progress  on  this  project. 

2.  On  July  7,  2010,  the  PI  visited  AFRL-Hanscom  and  discussed  with  Drs.  Shore  and  Yaghjian 
on  our  recent  progress  and  future  directions  for  this  project. 

3.  The  PI  attended  2011  Electromagnetics  Contractors  Meeting  organized  by  Dr.  Arje  Nachman 
from  AFOSR  between  January  4  and  5  and  met  Dr.  Yaghjian  and  other  researchers  in  this 
area. 

What  is  Next? 

1.  Design  antennas  using  DNG  metamaterials:  design  patch  antennas,  aperture  antennas  with 
excitation  of  leaky  waves,  optimize  microstrip  antenna  performance  (gain  and  bandwidth)  for 
different  dimensions,  thickness,  loss,  permittivity,  and  permeability. 

2.  Design  antennas  with  real  structures:  use  periodic  structures  for  substrate  and  superstrate. 
and  work  with  experimental  experts  to  make  such  antennas. 

3.  Reduce  loss  in  metamaterials  and  increase  antenna  efficiency. 

These  antennas  may  have  significant  impacts  on  applications  for  Air  Force. 
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